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ABSTRACT Black or dark brown (phaeoid) fungi cause cutaneous, subcutaneous, and systemic infections in KEYWORDS 

humans. Black fungi thrive in stressful conditions such as intense light, high radiation, and very low pH. Wangiella human 
{Exophiala) dermatitidis is arguably the most studied phaeoid fungal pathogen of humans. Here, we report our pathogenic 
comparative analysis of the genome of W. dermatitidis and the transcriptional response to low pH stress. This fungi 

revealed that W. dermatitidis has lost the ability to synthesize alpha-glucan, a cell wall compound many path- comparative 
ogenic fungi use to evade the host immune system. In contrast, W. dermatitidis contains a similar profile of chitin genomics 

synthase genes as related fungi and strongly induces genes involved in cell wall synthesis in response to pH RNA-Seq 

stress. The large portfolio of transporters may provide W. dermatitidis with an enhanced ability to remove harmful cell wall 
products as well as to survive on diverse nutrient sources. The genome encodes three independent pathways for biosynthesis 

producing melanin, an ability linked to pathogenesis; these are active during pH stress, potentially to produce horizontal gene 
a barrier to accumulated oxidative damage that might occur under stress conditions. In addition, a full set of transfer 
fungal light-sensing genes is present, including as part of a carotenoid biosynthesis gene cluster. Finally, we 
identify a two-gene cluster involved in nucleotide sugar metabolism conserved with a subset of fungi and 
characterize a horizontal transfer event of this cluster between fungi and algal viruses. This work reveals how 
W. dermatitidis has adapted to stress and survives in diverse environments, including during human infections. 
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A diverse array of primary and secondary diseases in humans is caused 
by a plethora of fungi with dark brown to black (phaeoid) vegetative 
growth at all stages of development (Matsumoto et al 1994). Tradi- 
tionally, these mycoses mostly have been associated with chronic 
dermatotropic forms of cutaneous and subcutaneous phaeohyphomy- 
cosis, chromoblastomycosis, and black-grain eumycotic mycetoma 
(Matsumoto et al 2011). Recently, however, an increasing frequency 
of systemic forms, particularly among patients with phaeohyphomy- 
cosis, has been reported in immunocompromised individuals and 
immunocompetent individuals (Zeng et al 2007). In extreme cases, 
the systemic forms may manifest as neurotropic and cerebral disease 
with no clear evidence of initial site of infection (Chang et al 2000; 
Sutton et al 1998). Although the distributions of the known patho- 
genic phaeoid fungi were once believed limited to tropical and sub- 
tropical regions, today these species and their potentially pathogenic 
relatives are being detected worldwide in widely diverse ecological 
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settings. No matter what the origin, whether infected human tissue, 
partners in lichens, saprophytes in leaf liter or on other woody 
materials, surfaces in the interior of the Chernobyl nuclear reactor, 
the extreme environments at polar regions, or bare rock at the top of 
mountains, the unifying characteristic of all phaeoid fungi is the 
deposition of 1,8-dihydroxynaphthalene melanin into its cell walls 
at all stages of development (Dadachova et al 2007; Feng et al 2001; 
Wheeler et al 2008; Schnitzler et al 1999). 

Among phaeoid fungi in general and human pathogenic phaeoid 
fungi in particular, Wangiella dermatitidis (also referred to as Exophiala 
dermatitidis) is the best known species for a number of reasons. The 
two most important of these are that it is a paradigm for phaeo- 
hyphomycosis and it is the model of choice for the more than 100 
other phaeoid fungi known to cause human disease (Matsumoto et al 
1994; Szaniszlo 2002, 2006). The latter is particularly important be- 
cause although each phaeoid agent may cause relatively few infections, 
worldwide each year the phaeoid fungi collectively give rise to large 
numbers of cases that are notoriously troublesome, often require sur- 
gery, and frequently are associated with failed antifungal therapy and 
death (Matsumoto et al 1994, 2011; Zeng et al 2007; Chang et al 
2000; Chua et al 2001; Filizzola et al 2003). Thus, the clinical di- 
versity, the large number of agents, and the impracticality of studying 
each of the many phaeoid agents of disease in great depth dictate that 
a model approach is necessary to gain insights into the varied biology 
of the ever-expanding repertoire of phaeoid fungi associated with 
human disease. 

A recent example of the need for an overall better understanding 
of the phaeoid fungi is related to the discovery by the medical 
community of the virtual absence of important biological knowledge 
of the phaeoid mold Exserohilum rostratum (Centers for Disease 
Control and Prevention 2012). This usually benign and relatively un- 
studied fungus as of 5 November 2012 had caused more than 749 
cases of phaeohyphomycotic meningitis and at least 61 deaths in 19 
states in the United States (Andes and Casadevall 2013) when it was 
injected as a contaminant of pharmaceutically compounded steroid 
solutions administered for pain management. In contrast with the 
cutaneous and subcutaneous forms of phaeohyphomycosis, surgery 
is not an option for these patients; instead, they require long-term 
antifungal therapy that can be dangerous and often fails. Just 10 yr 
previously, W. dermatitidis itself was the cause of a similar but smaller 
outbreak of fungal meningitis attributable to the injection of contami- 
nated steroids (Centers for Disease Control and Prevention 2002). Ad- 
ditional parallels can be extrapolated from studies of the W. dermatitidis 
model; 1,8-DHN melanin, a known virulence factor of W. dermatitidis, 
was also identified in E. rostratum and could similarly impart resistance 
to host killing and antifungal drug therapy (Paolo et al 2006; Al-Abdely 
2009; Andes and Casadevall 2013). In addition, multiple chitin syn- 
thases most likely are present in E. rostratum, which are also known 
to be important in the W. dermatitidis model for resistance to attack 
during infections by host immune systems and to antifungal therapeutic 
drugs, as well as cell wall strengthening and cell wall remodeling 
(Abramczyk et al 2009). Additionally, the regulation of its development 
from a spore to a hypha is likely controlled by pathways similar to one 
or more of those already studied in W. dermatitidis (Ye and Szaniszlo 
2000; Wang and Szaniszlo 2007; Liu et al 2008). This last hypothesis is 
plausible because even though W. dermatitidis is usually identified as 
a black yeast species, it is, in the manner of E. rostratum, in fact a con- 
idiogenous phaeoid mold with the additional growth characteristic of 
being polymorphic. The experimental examination of the variable mor- 
phology of W. dermatitidis has led to the discovery of numerous reg- 
ulatory aspects of its well-defined blastic, apical, and isotropic growth 



modes that produce myriad black forms of yeast, hyphae, and sclerotic 
(muriform) growth, respectively (Szaniszlo 2002; Abramczyk et al 
2009). The results have confirmed that this one model produces most, 
if not all, of the many growth forms of most other phaeoid agents of 
disease and therefore should have developmental pathways in common 
with any and all of the developmentally less versatile species. 

W. dermatitidis has been characterized as having remarkable ther- 
motolerance, halotolerance, pH tolerance, and hydrotolerance (Zalar 
et al 2011). Both thermotolerance and hydrotolerance are docu- 
mented by the frequent isolation of W. dermatitidis from environ- 
ments including hot tubs, saunas, steam baths, humidifiers, and 
dishwashers (Nishimura and Miyaji 1982; Matos et al 2002; Zalar 
et al 2011). In addition, this species was recently demonstrated to 
display enhanced growth induced by ionizing radiation (Dadachova 
et al 2007; Robertson et al 2012). This combination of tolerances is 
rarely observed in a single fungus and additionally makes W. dermatitidis 
an excellent system for studies of the extremophilic nature of fungi. 
Finally, it is important to note that no other phaeoid human pathogenic 
fungus besides W. dermatitidis is yet molecularly tractable and amena- 
ble to molecular transformation, targeted gene disruption, molecular 
gene transformation, and site-specific gene expression experiments 
(Abramczyk et al 2009; Szaniszlo 2002). 

The genome sequencing of W. dermatitidis provides the first op- 
portunity to examine the large group of human-associated fungi of the 
order Chaetotheriales within the subclass Chaetotheriomycetidae of 
the Eurotiomycetes class of the Ascomycota. To expand the usefulness 
of the W. dermatitidis model, we describe here the first comparative 
genomic analysis of a phaeoid pathogen. We compare this genome to 
the most closely related sequenced fungi, including other pathogens, 
and find differences in gene content related to cell wall structure and 
secondary metabolite production. To characterize the genes involved 
in responding to low pH stress, we performed deep transcriptomic 
profiling and found that cell wall genes are significantly induced by 
this stress, suggesting that remodeling the cell wall is a major factor in 
the stress resistance in this organism. 

MATERIALS AND METHODS 
Growth conditions for W. dermatitidis 

Genomic DNA was prepared from the wild-type strain of W. dermatitidis 
8656 (ATCC 34100, Exophiala dermatitidis CBS 525.76) cultured in 
YPD (2% peptone, 1% Bacto Yeast extract, and 2% dextrose) broth, 
with shaking at 200 rpm at 25° for 24 hr. For RNA-Seq experiments, 
wild-type cells were cultured in defined minimal medium (3 g/L 
NaN0 3 , 1 g/L K 2 PH0 4 , 1 g/L MgS0 4 0.7 H 2 0, 0.5 g/L KC1, 0.003 g/L 
thiamine, 5.3 g/L NH 4 C1, 30 g/L glucose) adjusted to pH 6 and pH 2.5, 
respectively, with shaking at 25° for 24 hr. Total RNA was extracted 
with RiboPureYM- Yeast Kit (Invitrogen, Carlsbad, CA). For the light 
effect experiment, the wild- type and the wdpksl mutant with a disrupted 
polyketide synthase 1 gene were grown on either YPD agar or CDY 
(Czapek Dox plus 0.1% yeast extract) agar. 

Genome sequencing assembly and repetitive 
element identification 

For genome sequencing, we constructed two libraries, 180 base 
fragment and 2-kb to 3-kb jumps, and sequenced both on the 
Illumina HiSequation 2000 platform. The 101 -bp Illumina reads were 
assembled using ALLPATHS-LG (Gnerre et al 2010) (build 18 May 
2011, 09:01:19 R37280) with default parameters using approximately 
100-fold depth of fragment reads and 50-fold depth of jumping reads. 
The assembly included 10 scaffolds totaling 26,350,763 bases and 238 
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contigs totaling 26,346,203 bases. Mitochondrial sequence was assem- 
bled in a separate 26-kb contig. Consensus de novo repetitive elements 
in the assembly were identified using Repeat Scout (Price et al. 2005) 
version 1.0.5; additional copies of consensus elements were identified 
in the assembly using RepeatMasker version open-3.2.8 (Smit et al). 
Individual elements were classified based on BLAST to the NCBI 
nucleotide collection datasets (nr/nt) or to RepeatMasker reference 
sets (nucleotide and protein). 

RNA-Seq and annotation 

For RNA-Seq, strand- specific libraries were constructed from polyA- 
purified RNA using the dUTP second strand marking method 
(Parkhomchuk et al. 2009; Levin et al 2010) as previously described 
(Cuomo et al 2012); 101 -bp reads were generated on the Illumina 
HiSequation 2000 platform. RNA-Seq was assembled with the Inch- 
worm component of Trinity (Grabherr et al 201 1) and processed with 
PASA (Haas et al 2008) to generate a set of transcripts for gene 
prediction. An initial gene set was generated by using EvidenceModeler 
(Haas et al 2008) to select the best gene call for a given locus from 
multiple gene prediction programs and PASA RNA-Seq transcripts as 
previously described (Haas et al 2008, 2011). The RNA-Seq tran- 
scripts were compared to the EvidenceModeler gene set and missing 
genes were manually reviewed and edited when possible. Potential 
misannotated genes were also reviewed, including transcripts with 
>15 exons, introns >1000 bases, transcripts with Ns (uncertain 
bases), and CDS overlaps of >50 bases. Genes with similarity to 
repeats were identified by TPSI (http://transposonpsi.sourceforge. 
net/) or BLAST hits to known repeat proteins; high confidence 
repeats were removed from the gene set. 

Functional annotation and comparative analysis 

We compared the genome of W. dermatitidis to those of 10 other 
related fungi, including Aspergillus nidulans, A. niger, A. fumigatus, 
Coccidioides immitis, Trichophyton rubrum, Thielavia terrestris, Myce- 
liophthora thermophila, Neurospora crassa, Saccharomyces cerevisiae, 
and Schizosaccharomyces pombe. We first extracted protein sequences 
from all 1 1 genomes and generated functional annotations with Inter- 
ProScan (Zdobnov and Apweiler 2001), Blast2GO (Conesa et al 
2005), and BLAST against KEGG and COG. Transporter family clas- 
sifications were based on the top BLAST hit (expect less than le — 15) 
to TCDB (Saier et al 2009). Transmembrane regions were predicted 
using TMHMM (version 2.0) (Krogh et al 2001) and signal peptide 
sequences were predicted using SignalP 4.0 (Petersen et al 2011). 

To build ortholog clusters, all protein sequences from the 1 1 fungi 
were compared using an all-against-all BLAST (requiring le— 5), then 
clustered with the OrthoMCL2 pipeline (Li et al 2003) to identify 
ortholog groups using default parameters. To infer a phylogeny for the 
11 fungi, we first concatenated alignments of single- copy orthologs 
present in all 11 species; a phylogeny was estimated using RAxML 
(Stamatakis 2006) version 7.3.3 with 1000 bootstrap replicates and 
model PROTGAMMABLOSUM62. Syntenic regions were identified 
using DAGchainer (Haas et al 2004) to compare linkage of orthologs. 

To identify expansions of predicted functional domains, protein 
domain profiles were compiled for all 11 fungi compared using 
InterPro domain hits. Significant expansions of a particular InterPro 
domain in W. dermatitidis as compared to the total count in all other 
fungi were detected using the hypergeometric distribution using 
q-values (Storey and Tibshirani 2003) to correct for multiple testing 
(Table SI). Similarly, to identify transporter family expansions, trans- 
porter family counts in W. dermatitidis and the sum of the counts in 



other fungi were compared using Fisher exact test (which has 
a hypergeometric distribution), also using a q-value correction. The 
two largest transporter families [major facilitator superfamily (MFS) 
and amino acid/polyamine/organocation (APC)] were also classified 
based on conservation as follows: core clusters (present in each of four 
phylogenetic groups: W. dermatitidis; Aspergillus species; Onygenales 
(T. rubum and C. immitis); all other fungi analyzed); shared clusters 
that are present in at least two groups; and unique clusters that are 
only found within each group. Genes involved in functional pathways 
(Supporting Information, Table S2 and Table S3) were mapped by 
identifying orthologs of previously described gene sets for cell wall 
biosynthesis (De Groot et al 2009), the cell wall integrity pathway 
(Futagami and Goto 2012; Levin 2005), the UDP-GlcNAc pathway 
(Milewski et al 2006), the HOG pathway (Hohmann et al 2007; 
Furukawa et al 2005), the calcium/calcineurin signaling pathway 
(Rispail et al 2009), and the pH signaling pathway (Andersen et al 
2009). Secondary metabolite gene clusters were identified using 
SMURF (Khaldi et al 2010). Candidate gliotoxin biosynthetic 
genes in W. dermatitidis were identified based on ortholog mapping 
from the A. fumigatus gene cluster (Gardiner and Howlett 2005). 

Phylogenetic analysis 

For GT2 and chitin synthases, proteins containing the GT2 
domain were aligned with MUSCLE (Edgar 2004) and the phylogeny 
estimated with RAxML version 7.0.4 (Stamatakis 2006) with the 
PROTGAMMAWAGF model and 1000 bootstrap replicates. Phylog- 
enies of the UNDG and chitin-synthase-like genes were inferred with 
RAxML version 7.3.3 with model PROTGAMMAWAG and 1000 
bootstrap replicates. For UNGD/UGD classification, a maximum likeli- 
hood tree was inferred using MEGA5 (Tamura et al 2011), with 1000 
bootstrap replicates for UDP-N-acetylglucosamine 6-dehydrogenases 
(UNGD) UDP-glucose 6-dehydrogenases (UGD), and several experi- 
mentally characterized bacterial nucleotide sugar dehydrogenases. Previ- 
ously classified bacterial sequences from the Integrated relational Enzyme 
database (http://ebi.ac.uk/intenz) (Fleischmann et al 2004) were com- 
pared to fungal genes; these include one UDP-N-acetylglucosamine 
6-dehydrogenase (UNGD; EC:1.1.1.136) from Pseudomonas 
aeruginosa, one UDP-N-acetyl-D-mannosamine dehydrogenase (UNMD; 
EC: 1.1. 1.336) from Escherichia coli, one UNMD from Methanococcus 
maripaludiSy one GDP-mannose 6-dehydrogenase (GMD; EC: 1.1. 1.132) 
from Pseudomonas aeruginosa, and two UDP-glucose 6-dehydrogenases 
(UGD; EC: 1.1. 1.22) from Pseudomonas aeruginosa. Phylogenies 
of the UDP-glucose-6-dehydrogenase and GT1 genes were inferred 
using maximum likelihood with MEGA5 with 1000 bootstrap 
replicates. 

Differential gene expression using RNA-Seq 

We used differential expression analysis scripts in the Trinity pipeline 
(Grabherr et al. 2011; Haas et al. 2013) to process RNA-Seq data 
generated from two pH conditions (pH 2.5 and pH 6), with three 
biological replicates from each condition. Briefly, we first extracted 
protein coding gene sequences from the W. dermatitidis genome 
sequence based on coordinates of gene models and added 100 bases 
of flanking sequence on each side to approximate UTRs. Then, the 
RNA-Seq reads from each of the six samples were aligned to the 
extracted coding sequences using bowtie (Langmead et al. 2009). 
The alignment files were used to quantify transcript abundances by 
RSEM (Li and Dewey 2011). Differential gene expression analysis 
was conducted using edgeR with TMM normalization (Robinson 
et al. 2009; Kadota et al. 2012) using a corrected p- value (Benjamini 
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and Hochberg 1995) cut-off of le— 10. Analysis of the data using MA 
plots did not show evidence of any apparent bias. To detect gene 
function enrichment under each of the two pH conditions, gene sets 
classified based on InterPro, COG, KEGG, and TCDB functional 
annotations were tested for enrichment using GSEA (Subramanian 
et al. 2005), calculating q-values to correct for multiple testing. We 
also tested manually curated gene sets responsible for cell wall bio- 
synthesis (Table S2) or belonging to pathways involved in cell stress 
responses (Table S3). 

RESULTS 

Generation of a highly complete genome and gene set 

We generated a high-quality genome assembly that comprises a small 
number of scaffolds. We sequenced the genome of W. dermatitidis 
using Alumina technology and assembled with ALLPATHS. The as- 
sembly consisted of 10 scaffolds with N50 of 3.62 Mb, totaling 26.4 
Mb, with only 0.02% of gap bases in scaffolds. Some scaffolds likely 
corresponded to complete chromosomes; the smallest two scaffolds 
were 7 kb and 255 kb, suggesting there are eight or fewer chromo- 
somes ranging from approximately 1.2 to 4.2 Mb. In support of this, 
we identified a candidate telomeric repeat array at 10 of the scaffold 
ends; three of the scaffolds (1.2, 1.5, and 1.6) contained repeat arrays 
at both ends, as would be expected for complete chromosomes. A 
candidate telomeric repeat unit of "TTTAGGG" was identified that 
contained one additional base compared to the "TTAGGG" repeat 
found in telomeres of many eukaryotes including fungi. This new 
7-mer repeat was present in 6 to 13 copies at the identified scaffold 
ends. Utilizing this highly complete assembly, a total of 9269 genes 
were predicted using a combination of gene prediction programs and 
strand-specific deep sequencing of RNA (RNA-Seq) {Materials and 
Methods). A total of 8826 genes (95.2%) were supported by RNA- 
Seq (FPKM >0.5). 

W. dermatitidis shares the largest amount of syntenic gene blocks 
with other ascomycetes of the Eurotiomycetes class. An average of 
3.4 Mb of W. dermatitidis sequence was found in syntenic regions 
with each of the three aspergilli examined (Materials and Methods), 
with A. fumigatus sharing the highest amount of syntenic regions and 
genes. A slightly lower total of 2.9 Mb of sequence was found in 
syntenic regions with the two Onygenales species examined; approx- 
imately half (53%) of the 612 genes in these syntenic regions were also 
found in syntenic regions shared by W. dermatitidis and the 
Eurotiales. 

The assembly contained a small fraction of repetitive sequence, 
identified de novo based on sequence similarity to other regions of the 
genome (Materials and Methods). A total of 86 repeat consensus 
elements covered 2.9% of the genome. More than half of these sequen- 
ces (64%) shared similarity with characterized transposable elements. 
Most of the identified transposable elements in W. dermatitidis were 
retroelements (84%), with a 3-to-l ratio of Gypsy to Copia types. 
These are the most common types of LTR retroelements found in 
fungal genomes, and typically Gypsy elements have higher copy num- 
bers (Muszewska et al. 2011). 

Conservation of genes involved in mating and meiosis provides 
evidence for an as-yet-undetected sexual cycle (Schurko and Logsdon 
2008). The MAT locus in the sequenced isolate of W. dermatitidis 
contains a HMG gene (HMPREF1120_08862) flanked by SLA2 
(HMPREF1120_08859) and APN2 (HMPREF1120_08863), as in 
other Pezizomycotina ascomycetes (Lee et al. 2010). In addition, many 
of the genes defined as a core meiotic set are conserved in the genome; 
of 12 core meiotic genes, 9 are present in W. dermatitidis. The three 



missing genes, DMC1, HOP2 y and MND1, are also missing from 
organisms with known sexual cycles, including Caenorhabditis elegans 
and Drosophila melanogaster (Schurko and Logsdon 2008). These 
three genes are present in all the Eurotiales and Onygenales compared 
here, suggesting a relatively recent loss of these three genes in the 
W. dermatitidis lineage. 

Gene conservation and gene family evolution 

We compared gene conservation between W. dermatitidis and related 
fungi to characterize the unique content of W. dermatitidis. Within 
the Eurotiomycetes, W. dermatitidis is the first sequenced repre- 
sentative of the black yeasts that is part of the subclass Chaeto- 
thyriomycetidae, whereas aspergilli and Onygenales are part of the 
Eurotiomycetidae. We compared W. dermatitidis to three aspergilli 
(Aspergillus nidulans, A. niger, and A. fumigatus) and two Onygenales 
(Coccidioides immitis, Trichophyton rubrum), as well as outgroup 
fungi representing thermophiles (Thielavia terrestris, Myceliophthora 
thermophila) and three well-studied ascomycete models (Neurospora 
crassa, Saccharomyces cerevisiae, Schizosaccharomyces pombe) (Figure 1). 
A total of 2626 W. dermatitidis genes were conserved across all 11 
fungi compared, and a larger set of 4755 genes was conserved between 
W. dermatitidis and at least one other of the species examined. Despite 
the high divergence from the other fungi compared here, W. dermatitidis 
contained a similar number (1888) of species-specific genes as found 
in other fungi (Figure 1). These include 58 paralogous clusters (cor- 
responding to 129 genes) specific to W. dermatitidis. Examining the 
species-specific genes, the most abundant protein domains identified 
belong to families of transporters (MFS and amino acid transporters), 
transcription factors, and short-chain dehydrogenases. 

To examine differences in predicted functional content, we 
identified protein domains enriched in W. dermatitidis compared to 
other ascomycete fungi. A total of 1 1 InterPro domains were enriched 
in the genome of W. dermatitidis (Figure 2). These included major 
facilitator type (MFS) transporters, saccharide dehydrogenases, kinases, 
transcription factors, enoylreductases (found in polyketide synthases 
and dehydrogenases), and acyltransferases. No InterPro domain fam- 
ilies are significantly depleted in W. dermatitidis compared to those in 
the other fungi examined. 

To further examine transporter family expansions, we annotated 
transporter subfamilies in all 11 compared fungi using the curated 
transporter database TCDB (Saier et al. 2009). This analysis identified 
several transporter families that are significantly expanded (Fisher 
exact test, q-value < 0.05) in W. dermatitidis; these include the 
large MFS and amino APC transporter families, as well as sodium- 
dependent solute transporters, fatty acid transporters, and major in- 
trinsic protein transporters (Figure 3). The MFS and APC transporter 
families are expanded in W. dermatitidis as well as in the three asper- 
gilli relative to the other compared fungi (Figure 2). To examine if the 
expansions in W. dermatitidis and the Aspergillus species are indepen- 
dent, we analyzed the conservation of these two largest expanded 
families in W. dermatitidis, MFS (237 members) and APC (52 mem- 
bers) (Materials and Methods) (Figure SI). Conserved transporter clus- 
ters found in all species groups do not display expansions in 
W. dermatitidis or the aspergilli. In contrast, MFS and APC clus- 
ters conserved in a subset of species show significant copy number 
increases in these species (Figure SI A); transporter families shared 
between W. dermatitidis and Aspergillus species account for most 
of higher total number. For example, W. dermatitidis has a total of 
41 MFS transporter clusters that are shared with one other fungal 
group, of which 33 are shared with only the aspergilli. In addition, 
both species contain a large number (50 in W. dermatitidis and 100 



564 I Z. Chen et al. 



:~-ZLG3' Genes | Genomes | Genetics 



0.1 



97 
100 



W. derma titidis 
-A. nidulans 



1 A. niger 

1 — A. fumigatus 
— C. immitis 
T. rubrum 



100 1 — M. 
I — 7.1 



100 1 — M. thermophila 
terrestris 
N. crassa 

S. cerevisiae 



S. pombe 



■ Core 
I Conserved 
Species-specific 



5000 



10000 



15000 



Figure 1 Phylogenetic relation- 
ship and gene conservation of 
W. dermatitidis and 10 com- 
pared fungi. Phylogeny (left 
panel) is based on concate- 
nated MUSCLE alignments of 
2013 single copy genes. Species 
phylogeny was inferred using 
RAxML (BLOSUM62 matrix and 
PROTGAMMA model) with 1000 
bootstrap replicates. Ortholog con- 
servation (right panel) highlights 
genes conserved in all species 
(core, blue), genes conserved in 
at least two species (conserved, 
red), and genes unique to a given 
species (species-specific, green). 



in the aspergilli) of species- specific transporter clusters. These 
comparisons suggest that both lineage-specific expansion and copies 
shared by only W. dermatitidis and the aspergilli duplications 
contribute to the MFS and APC transporter family expansions 
in these species. 

A total of 12 acyltransferase proteins are found in W. dermatitidis, 
nearly twice the number as in all other analyzed fungi (Figure 2). 
Whereas four of the acyltransferases are conserved between 
W. dermatitidis and other fungi, the other eight copies appear 
species-specific. These proteins have between three and eight 
predicted transmembrane regions, suggesting they are likely embedded 
in the membrane. In eukaryotes, the function of this domain is un- 
known, but the membrane-embedded OatA acyltransferase from the 
pathogenic bacterium Staphylococcus aureus is suggested to provide 



resistance to lysozymes induced during early host defense response 
(Bera et al 2006). 

To identify genes that may be involved in pathogenesis and other 
unique phenotypic properties of W. dermatitidis, we compared the Inter- 
Pro profile to other human-associated ascomycete fungi (Figure 2). 
Although found living on human skin, W. dermatitidis most commonly 
forms cysts or abscesses more deeply among living cells of cutaneous 
and subcutaneous tissues; T. rubrum, the cause of athlete's foot, utilizes 
mainly the keratinized dead cells of the outermost surface of skin as 
well as the keratin of nail. To determine if there are any gene changes 
common to these two different types of dermal associations, we identi- 
fied InterPro domains enriched in both W. dermatitidis and T. rubrum. 
Significant domain enrichments common to the two species include 
three classes of kinases: tyrosine; serine/threonine; and aminoglycoside 
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Figure 2 Protein domains enriched 
in W. dermatitidis in comparison 
to other fungi. Domain counts in 
W. dermatitidis were compared to 
the other 10 fungi and enriched 
domains were identified using 
the hypergeometric distribution. 
Domains are sorted by p-value 
and then grouped by category 
similarity. All domains shown 
are significantly enriched, with 
q-values <0.05. 
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Figure 3 Transporter family expansion and expression enrichment under stress conditions. (A) Heatmap showing expansion of transporter 
families in W. dermatitidis compared to other fungi. Families that are significantly expanded (q-value <0.05) are marked with a triangle on the left 
of the heatmap as follows: red, detected in comparison between W. dermatitidis compared to other fungi excluding Aspergillus species; green, 
W. dermatitidis and aspergilli compared to all other fungi; and blue, W. dermatitidis compared to aspergilli. (B) Enrichment of transporter families 
under pH (green bars) or radiation stress (blue bars). Significant enrichments noted as **q-value <0.05; ***q-value <0.01. 



phosphotransferases (Table SI). Major facilitator transporters and poly- 
ketide synthases are next most highly overrepresented. In addition, we 
found a depletion of domains involved in the degradation of plant 
components in W. dermatitidis (Table SI) similar to that found in 
Onygenales species (Sharpton et al. 2009; Desjardins et al. 2011; 
Martinez et al. 2012). Pectate lyase and cellulose binding domains 
are reduced in both fungi, although higher in number in W. dermatitidis 
than in T. rubrum (Table SI). The glycosyl hydrolase family 28 (GH28) 
involved in degrading pectin was present in 11 to 20 copies in our 
comparison group of aspergilli but was absent in W. dermatitidis 
as well as in T. rubrum and C. immitis (Table SI). The GH7 family 
involved in cellulose degradation was also absent from all three 
species but present in two to six copies in the other examined fungi 
(Table SI). 

Unlike T. rubrum and C. immitis, W. dermatitidis does not contain 
a significant enrichment of any protease family genes (Table SI). 
Whereas GH family genes as a whole are less abundant in T. rubrum 
and C. immitis, these species both contain a subtilisin S8 protease 
expansion absent in W. dermatitidis (Table SI). The expanded family 
of subtilisins in C. immitis, which may play a role in keratin degradation, 



are fewer in number in W. dermatitidis; the subtilisin protease 
domain (IPR000209, Peptidase S8/S53, subtilisin/kexin/sedolisin) 
is found in only three genes in W. dermatitidis compared to 21 
genes in C. immitis. Although an enhanced capacity to degrade 
proteins may compensate for the lack of saccharide degradation 
in some species (Desjardins et al 2011), the substrate preference in 
W. dermatitidis appears intermediate between protease and carbo- 
hydrate specialists. 

Potential for secondary metabolite production 

W. dermatitidis contains a subset of the secondary metabolite enzymes 
found in other fungi. There are only seven NRPS or NRPS-like pro- 
teins and four PKS or PKS-like proteins; no hybrid PKS-NRPS was 
identified. By comparison, aspergilli contain a larger portfolio of these 
enzymes, with an average number of proteins across eight species of 
26 PKS or PKS-like, 28 NRPS or NRPS-like, 2 hybrid PKS-NRPS, and 7 
DMATS (Khaldi et al 2010). Of the 11 W. dermatitidis enzymes, 7 are 
widely conserved in at least 8 of the other species in our comparative 
set. Of the four others, three are found in all three aspergilli and one 
(HMPREF1120_03318) is specific to W. dermatitidis. 
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Of the nine potential NRPS or PKS clusters containing accessory 
genes, one NRPS (HMPREF1120_02993) cluster shares similarity with 
the gliotoxin cluster of A. fumigatus. In W. dermatitidis, the gliP-like 
NRPS is conserved along with seven of the linked genes in the 
A. fumigatus cluster, although they are highly rearranged in gene 
order in W. dermatitidis (Figure S2). This cluster is not highly 
expressed or regulated by pH; in contrast, a different NRPS gene 
(HMPREF1120_04809) had 11 of the 19 clustered genes, including 
the NRPS, repressed at low pH (corrected p-value <le — 10). 
Gliotoxin is produced by species of Trichoderma and Penicillium, 
in addition to species of Aspergillus, but has not been reported in 
W. dermatitidis or other black yeasts. The presence and expression of 
this cluster in W. dermatitidis suggest it is capable of producing gliotoxin 
or another similar molecule; the phytopathogen Leptosphaeria maculans 
produces a structurally similar molecule, sirodesmin, from a homologous 
gene cluster (Gardiner and Howlett 2005). The immunosuppressive 
properties of gliotoxin may explain the decreased virulence of Asper- 
gillus gliP mutants (Sugui et al. 2007); the presence of this cluster in 
W. dermatitidis suggests it has the potential to play a similar role in 
this organism. 

Novel aspects of cell wall biosynthesis 

Analysis of cell wall genes revealed potentially unique features of the 
W. dermatitidis cell wall. Chitin synthase genes, which belong to the 
CAZy glycosyl transferase GT2 family, encompass several different 
subfamilies involved in chitin generation. We identified all seven pre- 
viously described chitin synthase genes in W. dermatitidis (Table S2). 
In addition to these seven WdCHS genes, we identified one novel 
chitin synthase-like gene (HMPREF1120_01791) with a significant 
(4e— 23) match to the chitin synthase domain (PF03142). Based on 
a phylogenetic analysis of CHS genes, this novel CHS-like gene is most 
closely related to CHS4 and CHS5 but clearly separate from each 
group (Figure S3). This novel CHS has five predicted transmembrane 
regions, similar to the other known CHS. Further analysis of the 
conservation of this novel CHS-like gene is described below. 

Although many of the genes involved in cell wall biosynthesis are 
conserved in W. dermatitidis, we found most of the genes involved in 
a-glucan metabolism are absent. Notably, we found that W. dermatitidis 
has a reduced set of proteins for 1,3-a-glucan synthesis and process- 
ing; the Ags family of 1,3-alpha glucan synthases and the Agn 
family of 1,3-alpha-glucanases are not present in the W. dermatitidis 
genome. In contrast, both the Ags and Agn families are present in 
multiple copies in Aspergillus species (Table S2). Another class of 
enzymes, alpha amylases, can perform transglycosylation or hydrolysis 
of glycosidic linkages; many fungi have intracellular and secreted 
alpha-amylases (Van Der Kaaij et al 2007). A conserved intracellular 
alpha-amylase is present in W. dermatitidis (HMPREF1120_03460) 
and most other fungi analyzed. The secreted amylases appear more 
diverse; in W. dermatitidis, this protein (HMPREF1120_08319) does 
not contain a secretion signal and is missing two domains (DUF1966 
and starch binding domain) that are conserved in other secreted alpha 
amylases. This pattern is also observed for T. rubrum; this species has 
also lost the Ags family of 1,3-alpha-glucan synthases, the Agn family 
of 1,3-alpha-glucanases, and the secreted alpha amylases. Only a single 
intracellular alpha-amylase is found in T. rubrum. 

Proteins involved in the regulation of chitin synthase activity, the 
export of chitin synthases, and chitin degradation (chitinase) are all 
conserved, as are proteins involved in glucan synthesis and processing 
(Table S2), one of which has been previously described (Szaniszlo 
et al. 1985). Two chitin deacetylases are conserved; however, neither 
is predicted to be secreted in W. dermatitidis. In other species these 



enzymes are secreted and act to modify chitin in the cell wall. In 
addition, W. dermatitidis lacks copies of two conserved chitosanases, 
similar to S. cerevisiae and S. pombe. This shift in the content of cell 
wall proteins suggests that 1,3-a-glucan synthesis and processing in 
W. dermatitidis is likely reduced compared to that in the aspergilli and 
other filamentous fungi and suggest that cell wall composition may 
also vary for other components. 

Examining the conservation of other genes involved in cell wall 
structure revealed expansions specific to W. dermatitidis. GH17 trans- 
glycosylases are highly conserved across the fungi we compared; in 
addition to the conserved copies, two additional paralogs are present 
in W. dermatitidis (Table S2). These extra copies could enable addi- 
tional cell wall cross-linking, which may strengthen the cell wall and 
increase survival during stressful conditions. The GT2 glycosyl trans- 
ferase family is also expanded in W. dermatitidis; 10 copies are pres- 
ent, which is significantly more than in the other fungi compared 
(hypergeometric test, q-value=0.11) and more than has been reported 
to date in any other genome. The aspergilli have the next highest 
count of seven GT2 genes. Phylogenetic and domain structure 
analysis revealed that the three GT2 genes (HMPREF1120_04699, 
HMPREF1120_5299, HMPREF1120_05696) contain a domain for 
cellulose synthase (IPR005150) (Figure S3). These proteins share se- 
quence similarity with homologs in the three aspergilli we examined, 
including the celA gene (AN8444) of A. nidulans, which was pre- 
viously suggested to play a role in (3 -glucan synthesis (De Groot et al. 
2009). 

Glucanases containing the GH16 module provide an example of 
domain shuffling and potential cell wall innovation. In W. dermatiti- 
dis, one GH16 gene (HMPREF1120_03145) also contains multiple 
LysM domains; the combination of these domains has not been pre- 
viously detected. The LysM domain has recently been linked to path- 
ogenesis in plant pathogens (De Jonge et al. 2010) and mycoparasites 
(Gruber et al. 2011), and it is highly expanded in dermatophytes 
(Martinez et al. 2012). In all other filamentous fungi, the LysM do- 
main is found in combination with a GH18 domain; however, such 
a protein is not present in W. dermatitidis. Although the combination 
of domains found here is not necessarily pathogenesis-related, it again 
indicates that cell wall structure in W. dermatitidis is likely different, 
because the LysM could attach the GH16 to the cell wall, creating 
different linkages. 

Other expanded domain classes include metabolic functions that 
may be involved in modifying the cell wall of W. dermatitidis. Glucose 
ribitol dehydrogenase, which breaks down sugar ends that lack alco- 
hols, is found in high numbers in W. dermatitidis and aspergilli. 
Ribitol is found as a precursor in capsule formation in bacteria and 
may protect some fungi against osmotic stress (Pfyffer and Rast 1980). 

Horizontal transfer of two-gene cluster between fungi 
and algal viruses 

The novel CHS-like gene (CHSL) and the adjacent gene 
(HMPREF1120_01790), a predicted UDP-N-acetyl-D-glucosamine 
6-dehydrogenase (UNGD) (Figure S4), form a two-gene cluster with 
an unusual conservation pattern. This gene pair may work in concert; 
UNGD family enzymes provide precursors for glycosyltransferase 
enzymes. Although this gene pair is not conserved in the other fungi 
we compared, both genes are frequently observed in species from the 
Sordariomycetes class, including Fusarium fujikuroi, F. oxysporum, 
and F. graminearum, two Aspergillus species (A. flavus and A. oryzae), 
and three Ustilagomycotina species within the Basidiomycota, includ- 
ing Ustilago hordei (Figure 4). These genes are also adjacent in these 
other fungal genomes, such as the gene pair FOXB_ 16531 and 
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FOXB_16532 from R oxysporum and UHOR_08674 and UHOR_08675 
from U. hordei. This gene pair is present in multiple copies at unlinked 
loci in some Sordariomycetes; phylogenetic analysis suggests that lineage 
duplication events resulted in these extra copies (Figure 4). Whereas this 
gene pair is found in a subset of Aspergillus species, including A. flavus 
and A. oryzae, these genes are not part of a larger syntenic region shared 
with W. dermatitidis. The genes flanking this cluster are conserved 
within, but not between, species groups (Figure 4). 

In addition to fungal species, this pair of linked genes is also found 
in several chlorella viruses (Figure 4), a group of large viruses infecting 
green algae (Jeanniard et al 2013). A more distantly related gene pair 
is also found in the lysogenic brown algal phaeovirus EsV-1 and in the 
host Ectocarpus siliculosus y with a suggested function of either an 
alginate or a chitin synthase (Delaroque et al 2001). These green 
and brown algal viruses belong to two distantly related genera of 
the family Phycodnaviridae. Phylogenetic analysis places both families 
of viral genes outside the fungal gene clade for UNGD; the CHSL 
proteins form two clades, each with a group of either brown or green 
algal viral sequences at the base. In the clade without W. dermatitidis, 



a third gene, a polysaccharide deacetylase, is part of the gene cluster in 
nearly all species including the green algal viral genomes (Figure 4C). 
This analysis suggests independent horizontal transfer of this gene 
cluster from two distantly related groups of viruses. 

UNGDs and the related subfamily UDP-glucose 6-dehydrogenase 
(UGD; EC: 1.1. 1.22) belong to the UDP-glucose/GDP-mannose de- 
hydrogenase family, which converts different UDP-sugar monomers to 
acidic forms. The UNDG class generates UDP-N-acetylglucosarnmuronic 
acid (UDP-GlcNAcA) and the UGD class generates UDP- glucuronic 
acid (UDP-GlcA). We identified a candidate UGD protein in W. 
dermatitidis (HMPREF1120_03278) based on shared domain 
structure (NAD binding domain, a central domain and a UDP 
binding domain) and phylogenetic analysis (Figure S4). This 
UGD in W. dermatitidis is adjacent to a glycosyl transferase group 
1 gene (GT1; HMPREF1120_03277). Searching the NCBI protein 
database with these two W. dermatitidis proteins identified con- 
served copies of this two- gene cluster in many Eurotiomycetes and 
Sordariomycetes species but in only a single Dothideomycete spe- 
cies. Phylogenetic trees are consistent between these paired UGDs 
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Figure 4 Phylogenetic relationship and linkage conservation of the UDP-N-acetyl-D-glucosamine 6-dehydrogenase (UNDG) (A) and chitin 
synthase-like (B) gene pair. Phylogenies (A and B) were estimated using RAxML (Stamatakis 2006) with model PROTGAMMAWAG and 1000 
bootstrap replicates. Taxonomic classification is noted by the colored boxes in (A) and (B). This gene pair is mainly found in Ascomycetes, 
predominantly Sordariomycetes (pink box), but is also observed in Eurotiomycetes (blue), Dothidiomycetes (brown). Outside Ascomycota copies 
are conserved in two Basidiomycete species (both Ustillaginomycetes, maroon) and outside fungi in brown algal viruses (brown) and green algal 
viruses (green). Homologs of the UNDG in the brown algal viruses are too divergent for phylogenetic analysis. (C) Genome organization of the 
gene pair. These genes are adjacent in all genomes, including brown algal viruses, although the relative transcriptional orientation varies (arrow 
direction). 
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and GT1 proteins (Figure S5), supporting a common origin of 
these two genes. The A. fumigatus ortholog, UGD1, is essential 
for the production of an acidic glycolipid not previously described 
for any fungal species (Fontaine et al. 2009). The conservation of 
orthologs of this gene in W. dermatitidis and other fungi suggests 
that they are also capable of producing acidic glycolipids. 

Conservation of iron acquisition pathways 

Multiple iron transport pathways are present in the W. dermatitidis, 
enabling acquisition of iron from the host during infection. Ferrox- 
idase and metalloreductase enzymes could provide reduced ferric iron 
to putative high affinity iron permeases. As in A. fumigatus (Haas 
2012), the genes encoding the ferroxidase FetC (HMPREF1120_01590 
and HMPREF1120_04510) and high-affinity iron permease FtrA 
(HMPREF1120_01589 and HMPREF1120_04509) are adjacent in 
the genome; notably, W. dermatitidis has two copies of this two- gene 
cluster. An alternative pathway for iron acquisition involves sidero- 
phores. The proteins involved in synthesis of the siderophore fusarine 
C (FsC) in A. fumigatus (Haas 2012) are all present in W. dermatitidis; 
the SidG protein capable of converting FsC to triacetylfusarine C is 
not conserved. For these siderophore genes, the genomic organization 
is only partially conserved between these species. Although the SidD and 
SidF genes are closely linked in both genomes (HMPREF1120_01440 
and HMPREF1120_001438, respectively), the SidA and SidC 
genes are adjacent in W. dermatitidis (HMPREF1120_07635 
and HMPREF1120_07636, respectively), whereas in A. fumigatus SidI 
and SidC are adjacent. 

Melanin biosynthesis pathways and regulation 

Melanin is a significant component of the phaeoid fungal cell wall 
because it helps defend against diverse environmental stresses such as 
UV, radiation, and oxidizing agents. In pathogenic fungi, melanin has 
been implicated in virulence as well as defense against the host immune 
system (Nosanchuk and Casadevall 2003; Jacobson 2000). Three differ- 
ent pathways contribute to melanin production: the DHN-melanin 
pathway; the DOPA-melanin pathway; and the L-tyrosine degradation 
pathway. 

The DHN-melanin biosynthetic pathway is the best characterized 
pathway for producing fungal melanin, a black pigment localized to 
cell walls (Wheeler et al 2008). In aspergilli, several enzymes in this 
pathway are developmentally regulated because melanin is produced 
in conidia but not typically in hyphae (Youngchim et al 2004), in- 
cluding a polyketide synthase (PKSP), scytalone dehydratase (ARP1), 
1,3,6,8-tetrahydroxynaphthalene reductase (ARP2), an alpha/beta hy- 
drolase (AYG1), and two multicopper oxidases (MCOs), a ferroxidase 
(ABR1), and a pigment MCO (ABR2) (Table 1). Homologs for all 
these proteins are identified in W. dermatitidis, including the known 
polyketide synthase (HMPREF1120_03173, WdPkslp) and the alpha/ 
beta hydrolase AYG1 (HMPREF1120_02312, WdYGl), both of which 
are required for melanin production in W. dermatitidis (Feng et al 
2001; Wheeler et al 2008). A second candidate alpha-beta hydrolase 
similar to AYG1 was also identified (HMPREF_00377, WdYG2); this 
gene is more closely related to the AYG1 orthologs in aspergilli. This 
suggests the functional copy of AYG1 in this pathway has been re- 
cently replaced in W. dermatitidis. The WdPKSl gene is conserved in 
other Eurotiomycetes and in N. crassa; however, these represent PKSs 
involved in production of other secondary metabolites. Both Arpl and 
Arp2 are present as a single ortholog in W. dermatitidis. In contrast, 
oxidase Abr2 and oxidase Abrl each has multiple homologous genes 
in the genome (Table 1). 



Despite the notable production of melanin by W. dermatitidis, the 
DHN-melanin pathway genes are not clustered in the genome (Table 1). 
The genes in this pathway are clustered in many other fungi, including 
A. fumigatus, A. clavatus, Penicillium marnejfei, and P. chrysogenum, 
with some rearrangements in the cluster; however, the genes are not 
clustered in two aspergilli, A. nidulans and A. niger (Woo et al. 2010). 
Despite the lack of clustering, high transcript levels were detected for all 
genes in the pathway. Genes for each step of this pathway are highly 
transcribed under both pH conditions, with average expression of the 
entire pathway (FPKM 391.9 under pH 2.5 and FPKM 150.4 under 
pH 6) greater than the average expression of all genes (FPKM 
114.8). Compared to radiation stress, low pH stress caused more 
extensive induction of the DHN pathway: 5 out of the 12 genes 
were upregulated (corrected p-value <1E— 10) under pH 2.5; with 
elevated levels of radiation, only one gene is upregulated by the 
same cut-off (see Transcriptional response to pH shift section). 

MCOs are over-represented in W. dermatitidis, including ABR- 
like pigment-related oxidases (see above) and laccases (Table 1). 
W. dermatitidis contains a total of 13 MCO genes. By comparison, 
A. fumigatus and A. nidulans have seven and eight copies, respectively, 
whereas A. niger has 16 copies. These MCO genes are classified into 
three subfamilies: fungal pigment MCO; fungal ferroxidase; and asco- 
mycete laccase (Ramos et al. 2011). The expansion in W. dermatitidis 
involves recent paralogous duplications in all three groups (Table 1). 
Of the 13 MCOs in W. dermatitidis, 12 are predicted to be secreted. 
Eight are highly upregulated after a shift to low pH (corrected p-value 
<1E— 4), suggesting that most of these copies could be functional. 

In addition to DHN melanin, many fungi also produce black or 
brown melanins through the DOPA-melanin pathway (Langfelder 
et al. 2003) in which tyrosinases or laccases generate dopaquinone 
either through hydroxylation of L-tyrosine or through oxidation of 
DOPA. The dopaquinone can auto-oxidize and polymerize to form 
melanins. W. dermatitidis also uses the DOPA-melanin pathway for 
melanin production (Paolo et al. 2006) and contains high numbers of 
both laccases (described above) and tyrosinases. Of the four tyrosinase 
genes present in the genome, two are conserved with Aspergillus 
homologs; the other two copies appear unique to W. dermatitidis 
and are highly upregulated (more than 10-fold) under pH 2.5 stress 
(Table 1). Two laccases are also upregulated by more than four-fold; 
in combination with the high tyrosinase expression, this suggests in- 
creased production of DOPA melanin at low pH. Similar to the DHN 
pathway, the DOPA pathway genes are also highly expressed under 
both pH conditions, with average FPKMs of 276.6 under pH 2.5 and 
92.1 under pH 6. 

An alternative pathway to produce a third type of melanin, using 
the L-tyrosine degradation pathway to produce pyomelanin as 
a side-product (Keller et al. 2011; Schmaler-Ripcke et al. 2009), 
is also conserved in W. dermatitidis (Table 1). The two enzymes 
involved in this pathway are 4-hydroxyphenylpyruvate dioxygenase 
(hppD) and homogentisate dioxygenase (hmgA). The hppD converts 
4-hydroxyphenylpyruvate to homogentisate, which can be converted 
to pyomelanin through oxidation and polymerization. Homogentisate 
can alternatively be degraded to other compounds by HmgA and two 
downstream enzymes, maleylacetoacetate isomerase (MaiA) and 
fumarylacetoacetate hydrolase (FahA). In response to low pH, we 
found that only hppD was significantly upregulated (more than 
four-fold); all other genes in this pathway did not significantly vary 
in transcript levels, suggesting that under low pH W. dermatitidis 
might increase pyomelanin production by selectively activating hppD 
in the L-tyrosine degradation pathway. In contrast to the results from 
pH experiments, most of the genes in this pathway were significantly 
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downregulated under radiation stress; the hppD gene was repressed 
more than 20-fold in radiated cells (Table 1), suggesting W. dermatitidis 
was less capable of pyomelanin production under radiation stress. 

Light-sensing and pigment production 

Fungi have multiple light sensors that respond to different light 
wavelengths. These include a phytochrome that senses red light, 
the white-collar complex that senses blue light, a cryptochrome/ 
photolyase that responds to both blue and UVA light, and a fungal 
opsin that senses green light. Although these components are 
missing in Onygenales species including T. rubrum and C. immitis 
(Idnurm et al 2010), the W. dermatitidis genome contains ortho- 
logs for each class of light sensing proteins (Table 2). Under radi- 
ation treatment, the UVA light sensor cryptochrome CryA gene 
(HMPREF1120_04357) is most significantly upregulated, with 
transcript levels increasing more than three-fold (Table 2) (Robertson 
et al 2012). 

The velvet family of regulatory proteins (VeA, VelB, VelC, VosA) 
plays a key role in coordinating development and secondary 
metabolism in response to light. In A. nidulans, VeA was found to 
interact with the light sensor complex (FphA-LreB-LreA) (Bayram 
et al 2008). VeA can also form the heterotrimeric velvet complex 
with VelB and LaeA (VelB-VeA-LaeA), which regulates development 
and secondary metabolism. VelB and VosA can form a heterodimer 
(VosA- VelB) that represses asexual development during vegetative 
growth or in the absence of light. Orthologs for all of these proteins 
are found in W. dermatitidis (Table 2). W. dermatitidis has one extra 
copy of the VelC; these two proteins share the velvet factor domain 
(PF 11754), but only one copy shares orthologs with the other fungi 
compared in our analysis. 

In addition to these photosensory proteins, an opsin-like protein is 
conserved within a carotenoid biosynthetic four-gene cluster (Table 2 
and Figure 5). The carotenoid pathway genes are not always clustered; 
in some other species such as F. fujikuroi y all genes are adjacent 



(Thewes et al 2005) but in other species only subsets are clustered. 
In T. rubrum, one locus contains three genes, with an inversion event 
relative to W. dermatitidis and F. fujikuroi (Figure 5); the opsin 
gene is found on a different scaffold in the T. rubrum assembly. In 
W. dermatitidis, carotenoid pigments are involved in light recep- 
tion and also protect against oxidative damage (Geis and Szaniszlo 
1984) and ionizing radiation (Robertson et al 2012). The conser- 
vation of this cluster in W. dermatitidis is striking because the 
melanins produced by this fungus also serve similar roles. Exposure 
to light stimulates biosynthesis of melanin and carotenoid. Colonies of 
the wild-type of W. dermatitidis show darker pigmentation under 
light than those grown in the dark, and colonies of the melanin- 
defective mutant wdpksl show an increased pink color under light, 
indicating carotenoid production (Figure 5). The coordination of 
these responses would enable protection and response to light- 
related stimuli. 

Transcriptional response to pH shift 

Previous studies demonstrated that stress conditions induce morpho- 
logical changes in W. dermatitidis cells. For example, under pH 2.5, 
the yeast form grows isotropically into enlarged forms that may 
become multicellular with unique characteristics, such as being 
multinucleate, having thick cell walls, and exhibiting melanin 
overproduction (Figure 6) (Szaniszlo et al 1993). We analyzed gene 
expression under two pH conditions, pH 2.5 and pH 6, to examine 
how W. dermatitidis responds to acidic pH stress. For each pH con- 
dition, RNA was prepared for three biological replicates and expres- 
sion levels measured using RNA-Seq. The biological replicates have 
highly correlated expression values, with an average Pearson correla- 
tion coefficient between replicates of 0.98 for pH 2.5 and 0.99 for pH 
6. There is also high correlation between the two pH conditions, from 
0.90 to 0.94, suggesting many genes are similarly expressed in both 
conditions. By comparing the RNA-Seq from the two pH conditions 
using edgeR (Robinson et al 2009), a total of 1291 genes were 



Table 2 Light-sensing genes and carotenoid biosynthetic pathway genes 



pH a Radiation 3 



Gene 


N. crassa 


W. dermatitidis 


Log 2 Fold 
Change 


p-value 


Log 2 Fold 
Change 


p-value 


Description 


Phy-1 


NCU04834 


HMPREF1120_ 


07867 


-0.71 


6.70E- 


-06 


-0.22 


1.22E- 


01 


Phytochrome, response to red light 


Phy-2 


NCU05790 




















Wc-1 


NCU02356 


HMPREF1120_ 


_06318 


-0.16 


3.71E- 


-01 


0.16 


2.68E- 


-01 


White-collar complex, response 


Vvd 


NCU03967 


















to blue light 


Wc-2 


NCU00902 


HMPREF1120_ 


04504 


-0.06 


7.70E- 


-01 


0.42 


2.18E- 


-03 




Phr 


NCU08626 


HMPREF1120_ 


04357 


-0.18 


3.69E- 


-01 


1.62 


4.73E- 


29 


Cryptochrome/photolyase, response 






















to blue and UVA light 


Nop-1 


NCU 10055 


HMPREF1120_ 


00264 


0.84 


8.65E- 


-08 


0.11 


4.35E- 


-01 


Fungal opsin, green light-sensing 


LaeA 


NCU00646 


HMPREF1120_ 


06677 


-0.85 


1.06E- 


-06 


-1.09 


2.29E- 


-14 


Velvet complex (VelBA/eA/LaeA), 


VeA 




HMPREF1120_ 


03349 


0.25 


1.38E- 


-01 


-0.01 


9.76E- 


-01 


coordinates development 


VelB 


NCU02775 


HMPREF1120_ 


00072 


-0.96 


9.76E- 


-10 


0.18 


2.12E- 


-01 


and secondary metabolism 






















in response to light 


VelC 


NCU07553 


HMPREF1120_ 


05878 


0.20 


2.64E- 


-01 


0.09 


5.80E- 


-01 


Velvet family protein with 






HMPREF1120_ 


02378 


1.29 


3.52E- 


-08 


0.00 


1.00E+00 


homology to VeA 


VosA 


NCU05964 


HMPREF1120_ 


_06091 


0.85 


5.08E- 


-04 


0.74 


1.54E- 


-04 




CarO 


NCU01735 


HMPREF1120_ 


_02861 


0.67 


2.19E- 


-04 


0.53 


3.16E- 


-03 


Opsin-like protein 






HMPREF1120_ 


_06324 


-0.59 


1.49E- 


-01 


0.10 


7.02E- 


-01 


CarB 


NCU00552 


HMPREF1120_ 


_02862 


-0.94 


3.60E- 


-09 


2.35 


9.29E- 


-71 


Phytoene dehydrogenase 


CarRA 


NCU00585 


HMPREF1120_ 


_02863 


0.23 


1.88E- 


-01 


1.98 


7.98E- 


-53 


Carotene cyclase 


CarX 




HMPREF1120_ 


_02864 


0.73 


1.58E- 


-05 


0.44 


1.17E- 


-03 


Carotenoid oxygenase 



Log 2 fold change and corrected p-values calculated by edgeR. Significant gene expression changes (corrected p-values <1E— 10) from both pH and radiation 
experiments are in bold. 
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F. fujikuroi 
W. dermatitidis 
T. rubrum 



B 




5 Kb 



i carotenoid oxygenase 
i carotene cyclase 
i phytoene dehydrogenase 
i opsin-like protein 




Figure 5 Effect of light on pigment production for 
W. dermatitidis. (A) Carotenoid biosynthesis gene clus- 
ter. All four genes are part of a single cluster in W. 
dermatitidis and in F. fujikuroi. The T. rubrum cluster 
does not include the CarO opsin protein, which is found 
at a different genomic location. (B and C) Wild-type 
(right sectors) and melanin defective mutant wdpksIA 
(left sectors) were incubated at 25° in the absence (B) or 
presence (C) of light for 2 weeks and grown on nutrient- 
limited CDY agar. Melanin (brown or black pigment) 
shown in the wild-type strain and carotenoid (pink pig- 
ment) shown in wdpksIA are both induced under light. 



identified as significantly differentially expressed (corrected p-value 
<le— 10) (Materials and Methods). These include 832 upregulated 
and 459 downregulated genes at pH 2.5 relative to pH 6. 

Analysis of the differentially expressed genes suggests that the 
major responses to low pH are to strengthen the cell wall, to pump out 
cations, and to slow metabolism (Table S3 and Table S4). To identify 
functional enrichment under low pH stress, we ran gene set enrich- 
ment analysis (GSEA) (Subramanian et al. 2005) using gene sets de- 
fined by protein domains, KEGG pathways, transporter families, and 
other custom sets (Materials and Methods) (Table S3 and Table S4). 
Genes induced under low pH stress include several involved in cell 
wall functions, including cell wall biosynthesis (including chitin and 
glucan processing and melanin biosynthesis) but not chitin degrada- 
tion (Figure 6C and Figure S6). The cell wall integrity pathway, in- 
cluding the transmembrane sensor and sensor-transducer proteins, is 
highly upegulated; in contrast, none of the genes in the high osmo- 
larity glycerol signaling pathway changes significantly in response to 
lowered pH (Table S4). Transporters (both MFS and APC) are also 
induced under low pH stress; the anion-cation symporters are most 
highly induced. Cytochrome p450 proteins, a diverse set of oxidizing 
enzymes, are also highly induced by pH and were previously noted to be 
highly induced by radiation stress (Robertson et al 2012). Other signif- 
icantly induced sets include melanin pathway genes, transcription fac- 
tors, and aspartate proteases. For genes repressed under pH 2.5, enriched 
categories include housekeeping functions, such as mRNA export pro- 
teins. This may reflect slower growth of these cells at low pH levels. 

The upregulation of cell-wall-related genes helps explain previous 
findings that W. dermatitidis cells develop thicker cell walls in low pH 
(Szaniszlo 2002, 2006; Szaniszlo et al 1978). Most of the genes in- 



volved in chitin synthesis are significantly upregulated under low pH 
stress. Five out of the seven known chitin synthases are induced 2.5- 
fold to 3.6-fold under pH 2.5 (Table S2). In contrast, genes involved in 
glucan synthesis are not as significantly upregulated as chitin synthesis 
genes. In addition, genes involved in glucan modification are signifi- 
cantly upregulated under pH 2.5. For example, a glycosyl hydrolase 17 
gene (HMPREF1120_08078) homologous to the S. cerevisiae putative 
glucanase SCW11 is upregulated approximately 100-fold at pH 2.5. 
This supports previous results that showed that under low pH stress, 
cells produce thicker cell walls by upregulating cell wall biosynthesis 
genes, and that cells are enriched with chitin and reduced in glucans 
under low pH stress (Szaniszlo et al 1983). 

Calcium concentration plays a critical role in regulating cellular 
growth and polymorphism, and transporters appear induced at low 
pH. Both acidic stress and depletion of Ca 2+ at neutral pH induce 
a transition of yeast form cells to isotropic or multicellular forms 
(Karuppayil and Szaniszlo 1997). In accordance with the phenotypic 
data, three genes encoding calcium transporters (HMPREF1120_08350, 
HMPREF1120_05235, and HMPREF1120_04272) are significantly 
upregulated under pH 2.5. These proteins belong to three different 
superfamilies: voltage-gated ion channels, Ca 2+ /Cation antiporter, and 
P-type ATPase, respectively. The increased transcription suggests that 
under acidic growth condition, cells actively regulate Ca 2+ in response 
to low pH stress using these calcium transporters and also confirms at 
the molecular level that calcium is a critical factor for controlling 
cellular morphogenesis in W. dermatitidis. 

At pH 2.5, cells grow isotropically and eventually produce 
multicellular forms (containing multiple nuclei and septa in one cell). 
The APSES transcription factor is a regulator of yeast-hyphal 
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Figure 6 Response of W. dermatitidis to low pH. (A and B) Light microscopy of W. dermatitidis at pH 2.5 (top panel, magnification 1200x) and pH 
6.5 (lower panel, magnification 1600x). (C) Functional category enrichment measured by GSEA. The x-axis of the MA-plot measures the log fold 
difference, and the y-axis shows log fold change between the two pH conditions. Four categories (chitin synthase, glucan processing, melanin, 
and cell wall integrity) of cell wall-related genes are significantly induced at low pH; one displayed category (chitin degradation) is not significantly 
induced. 



transitions and is required for hyphal formation in W. dermatitidis 
(Wang and Szaniszlo 2007). This gene (WdSTUA, HMPREF1120_05541) 
is conserved in A. nidulans (stuA y ANID_05836), N. crassa (Asm-l y 
NCU01414), S. cerevisiae (PHD1,YKL043W; SOK2,YMR016C), as well 
as three other genomes in our data set (T. rubrum, C. intimitis, and A. 
fumigatus). At low pH, the transcript level of WdSTUA is significantly 
induced (corrected p-value=le— 7) by nearly two-fold. Fungal morpho- 
genesis also requires motor proteins, including kinesins and myosins; 
seven kinesins, one dynein, and four myosins were found in the 
W. dermatitidis genome. The myosin 17 family includes two chitin syn- 
thases (Chs5p and Chs7p) with myosin head motor domains. Among 
these motor proteins, the kinesin-7 (KipA, HMPREF1120_04488) and 
the two Chs proteins with motor domains (Chs5, HMPREF1120_08776; 
Chs7, HMPREF1120_08777) are significantly upregulated (corrected p- 
value <4e— 17) under low pH. Under radiation, when cell morphology 
does not shift to isotropic growth, most of these genes are significantly 
downregulated and only Kinesin-3 (UncA, HMPREF1120_07465) is 
significantly upregulated (corrected p-value=4.5e— 9). 

DISCUSSION 

This analysis of the genome and transcriptional response to pH of 
W. dermatitidis is enabled by a high-quality assembly and gene set. 
Most of the eight large assembled scaffolds likely correspond to 



complete chromosomes; pulsed field gel analysis has suggested 
a karyotype of four chromosome bands, of which two appear to 
be doublets (Zheng et al 2006). Supporting this, the sequence of 
just more than half of the ends of the eight large scaffolds consists 
of a newly identified telomeric repeat array. The three scaffolds 
with repeat arrays at both ends are likely complete chromosomes. 
In addition, there are a few internal gaps in this assembly. Gaps 
between contigs occupy only 0.02% of the total scaffold length. The 
deep-coverage RNA-Seq provided evidence for a highly complete 
and accurate gene set. 

Although many other fungal genomes have been sequenced to date, 
this analysis represents the first detailed genomic analysis of a member 
of the Chaetotheriomycetidae subclass of fungi. Within the Euro- 
tiomycetes class, distantly related fungi include aspergilli and Onyge- 
nales species, although W. dermatitidis and other related black 
fungi show a closer relationship to the fungi in some lichens (James 
et al. 2006). The sequencing of additional genomes within the 
Chaetotheriomycetidae and related lichen species will enable more 
detailed comparisons of how this group evolved in very different 
ecological niches. 

The transcriptional response to pH has also been investigated for 
other fungi, including A. niger (Andersen et al 2009) and C. albicans 
(Bensen et al 2004). These studies have suggested that the pal/pacC 
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pH signaling pathway is induced under alkaline pH conditions. Con- 
sistent with this, we found that the pal/pacC pathway in W. dermatitidis 
is not induced under pH 2.5 (Table S3), suggesting that this pathway is 
most likely specific for more neutral and alkaline pH response and is 
conserved in many fungi. In contrast, pathways that respond to low pH 
stress have not been well-defined. In this study, we identified several cell 
wall stress response pathways in W. dermatitidis based on ortholog 
mapping, including the cell wall integrity pathway (CWI), HOG signal- 
ing pathway (HOG), and Ca 2+ /calcineurin signaling pathway. We also 
included the pathway that synthesizes UDP-GlcNAc (from fructose- 
6-phosphate), the substrate for chitin synthesis. Examining the expres- 
sion changes of these pathways under low pH stress highlights that the 
CWI pathway is the major stress response pathway induced by low 
pH stress; the entire pathway is induced, with only the exception of 
Rom2, which is slightly downregulated (Table S3). Seven out of the 12 
genes are upregulated by approximately two-fold or more (corrected 
p-values <1E— 10). The HOG pathway is also upregulated, but none 
of these genes has a significance level of corrected p-values <1E— 10. 
The UDP-GlcNAc pathway is also slightly upregulated, consistent 
with our finding that chitin synthesis activity is upregulated under 
low pH (Table S2). Six out of the seven known chitin synthases, as well 
as genes involved in regulating chitin synthase activity, are induced, and 
genes involved in chitin degradation are repressed (Table S2). This 
suggests that W. dermatitidis cells increase their chitin content under 
low pH stress. Genes involved in beta-glucan processing and modifi- 
cation are also upregulated, including the GAS family of beta- 1,3- 
glucanosyltransferases (HMPREF1120_07283, HMPREF1120_01763, 
and HMPREF1120_03477, orthologs of yeast GAS family) that are 
involved in cell wall modification and maintenance. 

Whereas the genes involved in the DHN-melanin pathway are 
present in W. dermatitidis, they are not clustered in the genome as 
they are in many other fungi. In some species, including A. fumigatus, 
P. marneffeiy and Alternaria brassicicola, the DHN-melanin pathway 
genes are clustered (Woo et al. 2010; Tsai et al. 1999; Kimura and 
Tsuge 1993). One possible explanation for the existence of dispensable 
metabolic pathway gene clusters is that the close linkage of pathway 
genes may benefit the regulation of pathway expression (Keller and 
Hohn 1997). However, break-up of gene clusters is not uncommon in 
fungi, which may reflect the lack of selection pressure required to 
maintain the cluster or to signal the first step in the eventual loss of 
the pathway (Keller and Hohn 1997). Constitutive production of the 
high amount of DHN-melanin throughout the life cycle of W. dermatitidis 
suggests that clustering may not be required for coordinated regula- 
tion of the genes in this pathway. However, upregulation of melanin 
biosynthesis genes and a number of cell wall genes under pH 2.5 
indicates coordination of different pathways. Thickening cell wall 
makes more space available for deposition of increased melanin. This 
synergic effect would allow melanized fungi to survive and to adapt 
better to extreme environments. 

In the context of the redundancy and importance of melanin for 
W. dermatitidis, the conservation of a complete set of light- sensing 
genes might seem paradoxical. Although surface melanin may shield 
some proteins from detecting light, the fact that light induces pig- 
ments such as DHN-melanin in its cell wall and carotenoid pigment 
in its membranes (Figure 5) suggests that light is able to penetrate 
these pigments and can still be sensed by the photoreceptor proteins. 
Alternatively, carotenoid pigments are more important at cell stages 
with lower melanin levels, such as during rapid growth at neutral pH. 
Therefore, overproduction of pigments in response to light is 
likely not shielding light but rather serving as scavengers to pro- 
tect against oxidative damages caused by light, or is involved in 



processes known to require light in other fungi, such as sexual 
reproduction, phototropism, and conidiation (Purschwitz et al 
2006). Because melanin and carotenoid are examples of secondary 
metabolites, it will be interesting to further investigate the effect of 
light on other secondary metabolism pathways, such as the one for 
the gliotoxin-like cluster identified in this genome. A previous 
study demonstrated that ionizing radiation also significantly 
upregulated gene expression of cryptochrome CryA, a well-known 
photoreceptor (Robertson et al. 2012). This suggests that this 
photosensory protein may also sense radiation. The identification 
of a complete set of fungal photosensory genes in W. dermatitidis 
suggests that, beyond the filamentous fungus N. crassa, this black 
yeast is another excellent model organism to study fungal photo- 
biology in addition to radiation biology. 

Nucleotide sugar dehydrogenases catalyze NAD/NADP-dependent 
oxidation of nucleotide-linked sugar monomers to produce acidic 
nucleotide sugars, and different subfamilies of these enzymes recognize 
different nucleotide sugars. We identified UDP -glucose dehydrogenase 
(UGD) and UDP-N-acetyl-D-glucosamine 6-dehydrogenase (UNGD) 
in W. dermatitidis genome; each gene is adjacent to a glycosyltrans- 
ferase gene (GT1 family and GT2 family, respectively). UGDs convert 
UDP-glucose (UDP-Glc) to UDP -glucuronic acid (UDP-GlcA) and 
UNGDs convert UDP-N-acetyl-D-glucosamine (UDP-GlcNAc) to 
UDP-N-acetyl-D-glucosaminuronic acid (UDP-GlcNAcA). In bacte- 
ria, UDP-GlcA and UDP-GlcNAcA are used for the synthesis of 
exopolysaccharides (EPSs) and lipopolysaccharides (LPSs) that are 
critical to bacterial virulence. The bacterial UNGD proteins are in- 
volved in the cell surface antigen biosynthesis process and are often 
located in gene clusters (Miller et al 2004; Zhang et al 2006) that also 
contain a GT family gene (Larkin and Imperiali 2009). In fungi, UGDs 
and UNGDs have not been well-characterized. The A. fumigatus 
UGD1 is the only one known so far that has been shown to be 
essential for the production of an acidic glycolipid (Fontaine et al 
2009). Conservation of the UGD-GT1 gene pair in many Eurotiomy- 
cetes and Sordariomycetes suggests these fungi can also produce acidic 
glycolipid in their cell membrane. UDP-GlcNAc is synthesized from 
fructose-6-phosphate through a four-gene pathway (Milewski et al 
2006), which is conserved in W. dermatitidis (Table S3). UDP- 
GlcNAc is the substrate for chitin synthases to make cell wall chitin; 
this substrate is also used for GPI-anchor biosynthesis (Kinoshita and 
Inoue 2000) and glycosylation of cell wall proteins (Jin 2012). How- 
ever, the role of the acidic form (UDP-GlcNAcA) in fungal cells is 
unclear at present. Under low pH stress, six of the seven known chitin 
synthases are upregulated; however, the UNGD and the CHSL genes 
are downregulated, with similar fold change under low pH and radi- 
ation stress (Table S2), suggesting that this gene pair may not function 
to produce cell wall inner layer chitin. Instead, they may function to 
modify cell surface structures during infection. 

The conservation of the UNGD-CSHL gene cluster in some fungi 
and algal viruses suggests horizontal transfer may have occurred 
between them. Previous work has shown that chlorella viruses have 
acquired and transferred genes involved in chitin metabolism to their 
algal hosts (Blanc et al 2010). Although their association with fungi 
has not been previously suggested, the conservation of the UNGD- 
CSHL gene pair suggests these viruses may interact directly with these 
fungal species. The presence of this gene pair in two groups of Phy- 
codnaviridae would have required independent transfer events from 
these two viral groups. 

Analysis of the W. dermatitidis genome suggests a shift in cell 
wall composition relative to related ascomycete fungi. The absence 
of enzymes involved in 1,3 alpha-glucan synthesis and processing 
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suggests this component is likely absent or highly reduced from the 
cell wall. This is in accord with a previous study examining cell wall 
polysaccharides, which noted the absence of alpha-glucan and chito- 
san from yeast and sclerotic bodies of W. dermatitidis (Geis 1981). 
This work also found that alpha- mannan formed the outer layer of 
both cell types, suggesting that these polysaccharides mediate surface 
interactions. 

Variation in cell wall composition predicted by the genome of 
W. dermatitidis likely affects interactions with the host immune 
system. In Histoplasma capsulatum and Magnaporthe oryzae, 
a-glucan helps protect cells from detection by the innate immune 
system (Rappleye et al. 2007; Fujikawa et al. 2012). In H. capsulatum, 
a-glucan may help mask P-glucans from detection by the mammalian 
receptor dectin-1. An alternative mechanism for masking P-glucans 
has been shown in A. fumigatus, in which hydrophobins help block 
immune cell recognition (Aimanianda et al. 2009). Chitin can also 
block host immune recognition in Candida albicans (Lee et al. 2012; 
Mora-Montes et al. 2011). Because W. dermatitidis appears to have 
lost the capacity to produce a-glucan, these cells may rely on chitin or 
a similar set of hydrophobin-like proteins to help protect cell wall 
P-glucans from detection or may alternatively produce more a-mannans. 
The upregulation of cell wall components at low pH suggests similar 
regulation could be adaptive for growing in a human host; increased 
chitin could help cells avoid detection by the host immune system. Al- 
ternatively, immune detection may be compromised, and the higher 
capacity of W. dermatitidis to evade stress may enable cells to survive. 
The ability of W. dermatitidis to cause severe infections in immunocom- 
petent hosts suggests it is able to overcome host defenses. 

ACKNOWLEDGMENTS 

We thank the Broad Institute Genomics Platform for generating all 
DNA and RNA sequence described here. The assembly and 
annotation were submitted to the NCBI nucleotide database under 
the accession number AFP AO 1000000. RNA-Seq differential expres- 
sion analysis of pH data were submitted to the NCBI GEO database 
(GSE51646). This work was supported by the National Human 
Genome Research Institute (grant number U54HG003067 to the 
Broad Institute), the Human Microbiome Project (grant U54HG004969), 
and Naval Research Laboratory core funds from the Office of Naval 
Research. P.J.S., Z.W., and C.A.C. designed the project. S.M.S. assembled 
genomes. S.G. and Q.Z. annotated genomes. Z.W. performed experi- 
ments. Z.C., D.A.M., Z.W., and C.A.C. analyzed data. Z.C., D.A.M., P.J.S., 
Z.W., and C.A.C. wrote the paper. 

LITERATURE CITED 

Abramczyk, D., C. Park, and P. J. Szaniszlo, 2009 Cytolocalization of the 
class V chitin synthase in the yeast, hyphal and sclerotic morphotypes of 
Wangiella (Exophiala) dermatitidis. Fungal genetics and biology. Fungal 
Genet. Biol. 46: 28-41. 

Aimanianda, V., J. Bayry, S. Bozza, O. Kniemeyer, K. Perruccio et al, 
2009 Surface hydrophobin prevents immune recognition of airborne 
fungal spores. Nature 460: 1117-1121. 

Al-Abdely, H. M., 2009 Phaeohyphomycosis: the dark question mark in 
clinical disease. J. Invasive Fungal Infect. 3: 82-88. 

Andersen, M. R., L. Lehmann, and J. Nielsen, 2009 Systemic analysis of the 
response of Aspergillus niger to ambient pH. Genome Biol. 10: R47. 

Andes, D., and A. Casadevall, 2013 Insights into fungal pathogenesis from 
the iatrogenic epidemic of Exserohilum rostratum fungal meningitis. 
Fungal genetics and biology. Fungal Genet. Biol. 61: 143-145. 

Bayram, O., S. Krappmann, M. Ni, J. W. Bok, K. Helmstaedt et al, 

2008 VelB/VeA/LaeA complex coordinates light signal with fungal de- 
velopment and secondary metabolism. Science 320: 1504-1506. 



Benjamini, Y., and Y. Hochberg, 1995 Controlling the false discovery rate: 
A practical and powerful approach to multiple testing. J. R. Stat. Soc. B. 
57: 289-300. 

Bensen, E. S., S. J. Martin, M. Li, J. Berman, and D. A. Davis, 

2004 Transcriptional profiling in Candida albicans reveals new adaptive 
responses to extracellular pH and functions for RimlOlp. Mol. Microbiol. 
54: 1335-1351. 

Bera, A., R. Biswas, S. Herbert, and F. Gotz, 2006 The presence of pepti- 
doglycan O-acetyltransferase in various staphylococcal species correlates 
with lysozyme resistance and pathogenicity. Infect. Immun. 74: 4598-4604. 

Blanc, G, G Duncan, I. Agarkova, M. Borodovsky, J. Gurnon et al, 
2010 The Chlorella variabilis NC64A genome reveals adaptation to 
photosymbiosis, coevolution with viruses, and cryptic sex. Plant Cell 22: 
2943-2955. 

Centers for Disease Control and Prevention, 2002 Exophiala infection from 
contaminated injectable steroids prepared by a compounding pharmacy — 
United States, July-November 2002. MMWR Morb. Mortal. Wkly. Rep. 51: 
1109-1112. 

Centers for Disease Control and Prevention, 2012 Multistate fungal men- 
ingitis outbreak — Interim guidance for treatment. MMWR Morb. 
Mortal. Wkly. Rep. 61: 842. 

Chang, C. L., D.-S. Kim, D. J. Park, H. J. Kim, C. H. Lee et al, 2000 Acute 
cerebral phaeohyphomycosis due to Wangiella dermatitidis accompanied 
by cerebrospinal fluid eosinophilia. J. Clin. Microbiol. 38: 1965-1966. 

Chua, J. D., S. M. Gordon, J. Banbury, G S. Hall, and G W. Procop, 2001 
Relapsing Exophiala jeanselmei phaeohyphomycosis in a lung-transplant 
patient. Transpl. Infect. Dis. 3: 235-238. 

Conesa, A., S. Gotz, J. M. Garcia-Gomez, J. Terol, M. Talon et al, 

2005 Blast2GO: a universal tool for annotation, visualization and 
analysis in functional genomics research. Bioinformatics 21: 3674- 
3676. 

Cuomo, C. A., C. A. Desjardins, M. A. Bakowski, J. Goldberg, A. T. Ma et al, 
2012 Microsporidian genome analysis reveals evolutionary strategies for 
obligate intracellular growth. Genome Res. 22: 2478-2488. 

Dadachova, E., R. A. Bryan, X. Huang, T. Moadel, A. D. Schweitzer et al, 
2007 Ionizing radiation changes the electronic properties of melanin 
and enhances the growth of melanized fungi. PLoS ONE 2: e457. 

Delaroque, N., D. G Muller, G Bothe, T. Pohl, R. Knippers et al, 2001 The 
complete DNA sequence of the Ectocarpus siliculosus virus EsV-1 ge- 
nome. Virology 287: 112-132. 

Desjardins, C. A., M. D. Champion, J. W. Holder, A. Muszewska, J. Goldberg 
et al, 2011 Comparative genomic analysis of human fungal pathogens 
causing paracoccidioidomycosis. PLoS Genet. 7: el002345. 

Edgar, R. C, 2004 MUSCLE: multiple sequence alignment with high ac- 
curacy and high throughput. Nucleic Acids Res. 32: 1792-1797. 

Feng, B., X. Wang, M. Hauser, S. Kaufmann, S. Jentsch et al, 

2001 Molecular cloning and characterization of WdPKSl, a gene in- 
volved in dihydroxynaphthalene melanin biosynthesis and virulence in 
Wangiella (Exophiala) dermatitidis. Infect. Immun. 69: 1781-1794. 

Filizzola, M. J., F. Martinez, and S. J. Rauf, 2003 Phaeohyphomycosis of the 
central nervous system in immunocompetent hosts: report of a case and 
review of the literature. Int. J. Infect. Dis. 7: 282-286. 

Fleischmann, A., M. Darsow, K. Degtyarenko, W. Fleischmann, S. Boyce 
et al, 2004 IntEnz, the integrated relational enzyme database. Nucleic 
Acids Res. 32: D434-D437. 

Fontaine, T., C. Lamarre, C. Simenel, K. Lambou, B. Coddeville et al, 
2009 Characterization of glucuronic acid containing glycolipid in As- 
pergillus fumigatus mycelium. Carbohydr. Res. 344: 1960-1967. 

Fujikawa, T., A. Sakaguchi, Y. Nishizawa, Y. Kouzai, E. Minami et al, 
2012 Surface a-l,3-glucan facilitates fungal stealth infection by inter- 
fering with innate immunity in plants. PLoS Pathog. 8: el002882. 

Furukawa, K., Y. Hoshi, T. Maeda, T. Nakajima, and K. Abe, 

2005 Aspergillus nidulans HOG pathway is activated only by two- 
component signalling pathway in response to osmotic stress. Mol. Microbiol. 
56: 1246-1261. 

Futagami, T., and M. Goto, 2012 Putative cell wall integrity sensor proteins 
in Aspergillus nidulans. Commun. Integr. Biol. 5: 206-208. 



G3'Genes | Genomes | Genetics 



Volume 4 April 2014 I Genomic and Transcriptomic Analysis of VV. dermatitidis I 575 



Gardiner, D. M., and B. J. Howlett, 2005 Bioinformatic and expression 

analysis of the putative gliotoxin biosynthetic gene cluster of Aspergillus 

fumigatus. FEMS Microbiol. Lett. 248: 241-248. 
Geis, P. A., 1981 Chemical Composition of the Yeast and Sclerotic Cell Walls of 

Wangiella dermatitidis. The University of Texas, Austin, TX. 
Geis, P. A., and P. J. Szaniszlo, 1984 Carotenoid pigments of the dematia- 

ceous fungus Wangiella dermatitidis. Mycologia 76: 268. 
Gnerre, S., I. Maccallum, D. Przybylski, F. J. Ribeiro, J. N. Burton et al, 

2010 High-quality draft assemblies of mammalian genomes from mas- 
sively parallel sequence data. Proc. Natl. Acad. Sci. USA. 108: 1513-1518. 

Grabherr, M. G., B. J. Haas, M. Yassour, J. Z. Levin, D. A. Thompson et al, 

2011 Full-length transcriptome assembly from RNA-Seq data without 
a reference genome. Nat. Biotechnol. 29: 644-652 

De Groot, P. W. J., B. W. Brandt, H. Horiuchi, A. F. J. Ram, C. G. de Koster 
et al, 2009 Comprehensive genomic analysis of cell wall genes in As- 
pergillus nidulans. Fungal genetics and biology. Fungal Genet. Biol. 46 
(Suppl 1): S72-S81. 

Gruber, S., G. Vaaje-Kolstad, F. Matarese, R. Lopez-Mondejar, C. P. Kubicek 
et al, 2011 Analysis of subgroup C of fungal chitinases containing 
chitin-binding and LysM modules in the mycoparasite Trichoderma 
atroviride. Glycobiology 21: 122-133. 

Haas, H., 2012 Iron — A key nexus in the virulence of Aspergillus fumigatus. 
Frontiers in microbiology 3: 28. 

Haas, B. J., A. L. Delcher, J. R. Wortman, and S. L. Salzberg, 

2004 DAGchainer: a tool for mining segmental genome duplications 
and synteny. Bioinformatics 20: 3643-3646. 

Haas, B. J., S. L. Salzberg, W. Zhu, M. Pertea, J. E. Allen et al, 

2008 Automated eukaryotic gene structure annotation using EViden- 
ceModeler and the Program to Assemble Spliced Alignments. Genome 
Biol. 9: R7. 

Haas, B. J., Q. Zeng, M. D. Pearson, C. A. Cuomo, and J. R. Wortman, 
2011 Approaches to fungal genome annotation. Mycology 2: 118-141. 

Haas, B. J., A. Papanicolaou, M. Yassour, M. Grabherr, P. D. Blood et al, 2013 De 
novo transcript sequence reconstruction from RNA-seq using the Trinity 
platform for reference generation and analysis. Nat. Protoc. 8: 1494-1512. 

Hohmann, S., M. Krantz, and B. Nordlander, 2007 Yeast osmoregulation. 
Methods Enzymol. 428: 29-45. 

Idnurm, A., S. Verma, and L. M. Corrochano, 2010 A glimpse into the basis 
of vision in the kingdom Mycota. Fungal Genet. Biol. 47: 881-892. 

Jacobson, E. S., 2000 Pathogenic roles for fungal melanins. Clin. Microbiol. 
Rev. 13: 708-717. 

James, T. Y., F. Kauff, C. L. Schoch, P. B. Matheny, V. Hofstetter et al, 
2006 Reconstructing the early evolution of Fungi using a six-gene 
phylogeny. Nature 443: 818-822. 

Jeanniard, A., D. D. Dunigan, }. R. Gurnon, I. V. Agarkova, M. Kang et al, 
2013 Towards defining the chloro viruses: a genomic journey through 
a genus of large DNA viruses. BMC Genomics 14: 158. 

Jin, C, 2012 Protein glycosylation in Aspergillus fumigatus is essential for 
cell wall synthesis and serves as a promising model of multicellular eu- 
karyotic development. Int. J. Microbiol. 2012: 654251. 

De Jonge, R., H. P. van Esse, A. Kombrink, T. Shinya, Y. Desaki et al, 

2010 Conserved fungal LysM effector Ecp6 prevents chitin-triggered 
immunity in plants. Science 329: 953-955. 

Kadota, K., T. Nishiyama, and K. Shimizu, 2012 A normalization strategy 

for comparing tag count data. Algorithms Mol. Biol. 7: 5. 
Karuppayil, S. M., and P. J. Szaniszlo, 1997 Importance of calcium to the 

regulation of polymorphism in Wangiella [Exophiala) dermatitidis. J.. 

Med.. Vet. Mycol.. 35: 379-388. 
Keller, N. P., and T. M. Hohn, 1997 Metabolic pathway gene clusters in 

filamentous fungi. Fungal Genet. Biol. 21: 17-29. 
Keller, S., J. Macheleidt, K. Scherlach, J. Schmaler-Ripcke, I. D. Jacobsen et al, 

2011 Pyo melanin formation in Aspergillus fumigatus requires HmgX 
and the transcriptional activator HmgR but is dispensable for virulence. 
PLoS ONE 6: e26604. 

Khaldi, N., F. T. Seifuddin, G. Turner, D. Haft, W. C. Nierman et al, 

2010 SMURF: Genomic mapping of fungal secondary metabolite clus- 
ters. Fungal Genet. Biol. 47: 736-741. 



Kimura, N., and T. Tsuge, 1993 Gene cluster involved in melanin biosyn- 
thesis of the filamentous fungus Alternaria alternata. J. Bacteriol. 175: 
4427-4435. 

Kinoshita, T., and N. Inoue, 2000 Dissecting and manipulating the pathway 
for glycosylphos-phatidylinositol-anchor biosynthesis. Curr. Opin. Chem. 
Biol. 4: 632-638. 

Krogh, A., B. Larsson, G. von Heijne, and E. L. Sonnhammer, 2001 Predicting 
transmembrane protein topology with a hidden Markov model: application 
to complete genomes. J. Mol. Biol. 305: 567-580. 

Langfelder, K., M. Streibel, B. Jahn, G. Haase, and A. A. Brakhage, 

2003 Biosynthesis of fungal melanins and their importance for human 
pathogenic fungi. Fungal Genet. Biol. 38: 143-158. 

Langmead, B., C. Trapnell, M. Pop, and S. L. Salzberg, 2009 Ultrafast and 
memory-efficient alignment of short DNA sequences to the human ge- 
nome. Genome Biol. 10: R25. 

Larkin, A., and B. Imperiali, 2009 Biosynthesis of UDP-GlcNAc(3NAc)A 
by WbpB, WbpE, and WbpD: enzymes in the Wbp pathway responsible 
for O-antigen assembly in Pseudomonas aeruginosa PAOl. Biochemistry 
48: 5446-5455. 

Lee, K. K., D. M. Maccallum, M. D. Jacobsen, L. A. Walker, F. C. Odds et al, 
2012 Elevated cell wall chitin in Candida albicans confers echinocandin 
resistance in vivo. Antimicrob. Agents Chemother. 56: 208-217. 

Lee, S. C, M. Ni, W. Li, C. Shertz, and J. Heitman, 2010 The evolution of 
sex: a perspective from the fungal kingdom. Microbiol. Mol. Biol. Rev. 74: 
298-340. 

Levin, D. E., 2005 Cell wall integrity signaling in Saccharomyces cerevisiae. 
Microbiology and molecular biology reviews. Microbiol. Mol. Biol. Rev. 
69: 262-291. 

Levin, J. Z., M. Yassour, X. Adiconis, C. Nusbaum, D. A. Thompson et al, 
2010 Comprehensive comparative analysis of strand- specific RNA se- 
quencing methods. Nat. Methods 7: 709-715. 

Li, B., and C. N. Dewey, 2011 RSEM: accurate transcript quantification 
from RNA-Seq data with or without a reference genome. BMC Bioin- 
formatics 12: 323. 

Li, L., C. J. Stoeckert, and D. S. Roos, 2003 OrthoMCL: identification of 
ortholog groups for eukaryotic genomes. Genome Res. 13: 2178-2189. 

Liu, H., D. Abramczyk, C. R. Cooper, Jr, L. Zheng, C. Park et al, 

2008 Molecular cloning and characterization of WdTUPl, a gene that 
encodes a potential transcriptional repressor important for yeast-hyphal 
transitions in Wangiella (Exophiala) dermatitidis. Fungal Genet. Biol. 45: 
646-656. 

Martinez, D. A., B. G. Oliver, Y. Graser, J. M. Goldberg, W. Li et al, 2012 
Comparative genome analysis of Trichophyton rubrum and related der- 
matophytes reveals candidate genes involved in infection.MmBio. 3: 
e00259-e00212. 

Matos, T., G. S. de Hoog, A. G. de Boer, I. de Crom, and G. Haase, 
2002 High prevalence of the neurotrope Exophiala dermatitidis and 
related oligotrophic black yeasts in sauna facilities. Mycoses 45: 
373-377. 

Matsumoto, T., L. Ajello, T. Matsuda, P. J. Szaniszlo, and T. J. Walsh, 1994 
Developments in hyalohyphomycosis and phaeohyphomycosis. J. Med. 
Vet. Mycol. 32 (Suppl 1): 329-349. 

Matsumoto, T., C. R. Cooper, and P. J. Szaniszlo, 2011 Chapter 84 - 
Chromoblastomycosis and phaeohyphomycosis, pp. 569-572 in Trop- 
ical Infectious Diseases: Principles, Pathogens and Practice, Ch 84, Ed. 3. 
W.B. Saunders, Edinburgh. 

Milewski, S., I. Gabriel, and J. Olchowy, 2006 Enzymes of UDP-GlcNAc 
biosynthesis in yeast. Yeast 23: 1-14. 

Miller, W. L., C. Q. Wenzel, C. Daniels, S. Larocque, J.-R Brisson et al, 

2004 Biochemical characterization of Wbp A, a UDP-N-acetyl-D-glucosamine 
6-dehydrogenase involved in O-antigen biosynthesis in Pseudomonas 
aeruginosa PAOl. J. Biol. Chem. 279: 37551-37558. 

Mora-Montes, H. M., M. G. Netea, G. Ferwerda, M. D. Lenardon, G. D. 
Brown et al, 2011 Recognition and blocking of innate immunity cells 
by Candida albicans chitin. Infect. Immun. 79: 1961-1970. 

Muszewska, A., M. Hoffman-Sommer, and M. Grynberg, 2011 LTR retro- 
transposons in fungi. PLoS ONE 6: e29425. 



576 I Z. Chen et a/. 



:~-ZLG3' Genes | Genomes | Genetics 



Nishimura, K., and M. Miyaji, 1982 Studies on a saprophyte of Exophiala 
dermatitidis isolated from a humidifier. Mycopathologia 77: 173-181. 

Nosanchuk, J. D., and A. Casadevall, 2003 The contribution of melanin to 
microbial pathogenesis. Cell. Microbiol. 5: 203-223. 

Paolo, W. F., Jr, E. Dadachova, P. Mandal, A. Casadevall, P. J. Szaniszlo et al, 

2006 Effects of disrupting the polyketide synthase gene WdPKSl in 
Wangiella [Exophiala] dermatitidis on melanin production and resistance 
to killing by antifungal compounds, enzymatic degradation, and extremes 
in temperature. BMC Microbiol. 6: 55. 

Parkhomchuk, D., T. Borodina, V. Amstislavskiy, M. Banaru, L. Hallen et al, 
2009 Transcriptome analysis by strand-specific sequencing of comple- 
mentary DNA. Nucleic Acids Res. 37: el23. 

Petersen, T. N., S. Brunak, G. von Heijne, and H. Nielsen, 201 1 SignalP 4.0: 
discriminating signal peptides from transmembrane regions. Nat. Meth- 
ods 8: 785-786. 

Pfyffer, G. E., and D. M. Rast, 1980 Systematic distribution of ribitol in the 
fungi. New Phytol. 85: 163-168. 

Price, A. L., N. C. Jones, and P. A. Pevzner, 2005 De novo identification of 
repeat families in large genomes. Bioinformatics 21(Suppl 1): i351-i358. 

Purschwitz, J., S. Miiller, C. Kastner, and R. Fischer, 2006 Seeing the rain- 
bow: light sensing in fungi. Curr. Opin. Microbiol. 9: 566-571. 

Ramos, J. A. T., S. Barends, R. M. D. Verhaert, and L. H. de Graaff, 

2011 The Aspergillus niger multicopper oxidase family: analysis and 
overexpression of laccase-like encoding genes. Microb. Cell Fact. 10: 78. 

Rappleye, C. A., L. G. Eissenberg, and W. E. Goldman, 2007 Histoplasma 
capsulatum a-(l,3)-glucan blocks innate immune recognition by the 
P-glucan receptor. Proc. Natl. Acad. Sci. USA 104: 1366-1370. 

Rispail, N., D. M. Soanes, C. Ant, R. Czajkowski, A. Grimier et al, 

2009 Comparative genomics of MAP kinase and calcium-calcineurin 
signalling components in plant and human pathogenic fungi. Fungal 
genetics and biology. Fungal Genet. Biol. 46: 287-298. 

Robertson, K. L., A. Mostaghim, C. A. Cuomo, C. M. Soto, N. Lebedev et al, 

2012 Adaptation of the black yeast Wangiella dermatitidis to ionizing 
radiation: molecular and cellular mechanisms. PLoS ONE 7: e48674. 

Robinson, M. D., D. J. McCarthy, and G. K. Smyth, 2009 edgeR: a Bio- 
conductor package for differential expression analysis of digital gene 
expression data. Bioinformatics 26: 139-140. 

Saier, M. H., Jr, M. R. Yen, K. Noto, D. G. Tamang, and C. Elkan, 2009 The 
transporter classification database: recent advances. Nucleic Acids Res. 
37: D274-D278. 

Schmaler-Ripcke, J., V. Sugareva, P. Gebhardt, R. Winkler, O. Kniemeyer 
et al, 2009 Production of pyomelanin, a second type of melanin, via the 
tyrosine degradation pathway in Aspergillus fumigatus. Appl. Environ. 
Microbiol. 75: 493-503. 

Schnitzler, N., H. Peltroche-Llacsahuanga, N. Bestier, J. Zundorf, R. Lutticken 
et al, 1999 Effect of melanin and carotenoids of Exophiala (Wangiella) 
dermatitidis on phagocytosis, oxidative burst, and killing by human 
neutrophils. Infect. Immun. 67: 94-101. 

Schurko, A. M., and J. M. Logsdon, 2008 Using a meiosis detection toolkit to 
investigate ancient asexual "scandals" and the evolution of sex. BioEssays 
30: 579-589. 

Sharpton, T. J., J. E. Stajich, S. D. Rounsley, M. J. Gardner, J. R. Wortman 
et al, 2009 Comparative genomic analyses of the human fungal 
pathogens Coccidioides and their relatives. Genome Res. 19: 1722-1731. 

Smit, A. F. A., R. Hubley, and P. Green RepeatMasker. Available at: http:// 
www.repeatmasker.org. version open-3.2.8. 

Stamatakis, A., 2006 RAxML-VI-HPC: maximum likelihood-based phylo- 
genetic analyses with thousands of taxa and mixed models. Bioinfor- 
matics 22: 2688-2690. 

Storey, J. D., and R. Tibshirani, 2003 Statistical significance for genomewide 
studies. Proc. Natl. Acad. Sci. USA 100: 9440-9445. 

Subramanian, A., P. Tamayo, V. K. Mootha, S. Mukherjee, B. L. Ebert et al, 
2005 Gene set enrichment analysis: a knowledge-based approach for 
interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 
102: 15545-15550. 

Sugui, J. A., J. Pardo, Y. C. Chang, K. A. Zarember, G. Nardone et al, 

2007 Gliotoxin is a virulence factor of Aspergillus fumigatus: gliP de- 



letion attenuates virulence in mice immunosuppressed with hydrocorti- 
sone. Eukaryot. Cell 6: 1562-1569. 
Sutton, D. A., M. Slifkin, R. Yakulis, and M. G. Rinaldi, 1998 U.S. case 
report of cerebral phaeohyphomycosis caused by Ramichloridium obo- 
voideum (R. mackenziei): Criteria for identification, therapy, and review 
of other known dematiaceous neurotropic taxa. J. Clin. Microbiol. 36: 
708-715. 

Szaniszlo, P. J., 2002 Molecular genetic studies of the model dematia- 
ceous pathogen Wangiella dermatitidis. Int. J. Med. Microbiol. 292: 
381-390. 

Szaniszlo, P. J., 2006 Virulence factors in black molds with emphasis on 
melanin, chitin, and Wangiella as a molecularly tractable model., pp. 407- 
428 in Molecular Principles of Fungal Pathogenesis. ASM Press, Wash- 
ington, D.C. 

Szaniszlo, P. J., P. Hsieh, and J. D. Marlowe, 1978 Induction and ultra- 
structure of the multiseptate (sclerotic) morphology in Phialophora der- 
matitidis. Mycologia 68: 117-130. 

Szaniszlo, P. J., M. S. Kang, and E. Cabib, 1985 Stimulation of P(l,3)glucan 
synthetase of various fungi by nucleoside triphosphates: generalized 
regulatory mechanism for cell wall biosynthesis. J. Bacteriol. 161: 
1188-1194. 

Szaniszlo, P. J., S. M. Karuppayil, L. Mendoza, and R. J. Rennard, 1993 Cell 
cycle regulation of polymorphism in Wangiella dermatitidis. Arch. Med. 
Res. 24: 251-261. 

Szaniszlo, P. J., P. A. Geis, C. W. Jacobs, C. R. Cooper, and J. L. Harris, 
1983 Cell wall changes associated with yeast to multicellular form 
conversion in Wangiella dermatitidis, pp. 239-244 in Microbiology, edited 
by D. Schlessinger. American Society of Microbiology, Washington, D.C. 

Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei et al, 

2011 MEGA5: molecular evolutionary genetics analysis using maxi- 
mum likelihood, evolutionary distance, and maximum parsimony 
methods. Mol. Biol. Evol. 28: 2731-2739. 

Thewes, S., A. Prado-Cabrero, M. M. Prado, B. Tudzynski, and J. Avalos, 
2005 Characterization of a gene in the car cluster of Fusarium fujikuroi 
that codes for a protein of the carotenoid oxygenase family. Mol. Genet. 
Genomics 274: 217-228. 

Tsai, H. F., M. H. Wheeler, Y. C. Chang, and K. J. Kwon-Chung, 1999 A 
developmentally regulated gene cluster involved in conidial pigment 
biosynthesis in Aspergillus fumigatus. J. Bacteriol. 181: 6469-6477. 

Van der Kaaij, R. M., S. Janecek, M. J. E. C. van der Maarel, and L. 
Dijkhuizen, 2007 Phylogenetic and biochemical characterization 
of a novel cluster of intracellular fungal alpha-amylase enzymes. 
Microbiology 153: 4003-4015. 

Wang, Q., and P. J. Szaniszlo, 2007 WdStuAp, an APSES transcription 
factor, is a regulator of yeast-hyphal transitions in Wangiella (Exophiala) 
dermatitidis. Eukaryot. Cell 6: 1595-1605. 

Wheeler, M. H., D. Abramczyk, L. S. Puckhaber, M. Naruse, Y. Ebizuka et al, 
2008 New biosynthetic step in the melanin pathway of Wangiella 
(Exophiala) dermatitidis: evidence for 2-acetyl-l,3,6,8-tetrahydroxynaphthalene 
as a novel precursor. Eukaryot. Cell 7: 1699-1711. 

Woo, P. C. Y., E. W. T. Tarn, K. T. K. Chong, J. J. Cai, E. T. K. Tung et al, 
2010 High diversity of polyketide synthase genes and the melanin 
biosynthesis gene cluster in Penicillium marneffei. FEBS J. 277: 3750- 
3758. 

Ye, X., and P. J. Szaniszlo, 2000 Expression of a constitutively active Cdc42 
homologue promotes development of sclerotic bodies but represses hy- 
phal growth in the zoopathogenic fungus Wangiella (Exophiala) derma- 
titidis. J. Bacteriol. 182: 4941-4950. 

Youngchim, S., R. Morris-Jones, R. J. Hay, and A. J. Hamilton, 

2004 Production of melanin by Aspergillus fumigatus. J. Med. Micro- 
biol. 53: 175-181. 

Zalar, P., M. Novak, G. S. de Hoog, and N. Gunde-Cimerman, 2011 Dishwashers- 
a man-made ecological niche accommodating human opportunistic fungal 
pathogens. Fungal Biol. 115: 997-1007. 

Zdobnov, E. M., and R. Apweiler, 2001 InterProScan-an integration plat- 
form for the signature- recognition methods in InterPro. Bioinformatics 
17: 847-848. 



G3'Genes | Genomes | Genetics 



Volume 4 April 2014 I Genomic and Transcriptomic Analysis of VV. dermatitidis I 577 



Zeng, J. S., D. A. Sutton, A. W. Fothergill, M. G. Rinaldi, M. J. Harrak et al, 
2007 Spectrum of clinically relevant Exophiala species in the United 
States. J. Clin. Microbiol. 45: 3713-3720. 

Zhang, H., Y. Zhou, H. Bao, and H. Liu, 2006 Vi antigen biosynthesis in 
Salmonella typhi: characterization of UDP-N-acetylglucosamine C-6 de- 
hydrogenase (TviB) and UDP-N-acetylglucosaminuronic acid C-4 
epimerase (TviC). Biochemistry 45: 8163-8173. 



Zheng, L., L. Mendoza, Z. Wang, H. Liu, C. Park et a/., 2006 WdChslp, 
a class II chitin synthase, is more responsible than WdChs2p (Class I) 
for normal yeast reproductive growth in the polymorphic, pathogenic 
fungus Wangiella (Exophiala) dermatitidis. Arch. Microbiol. 185: 
316-329. 

Communicating editor: A. Gasch 



578 I Z. Chen et a/. 



:~-ZLG3' Genes | Genomes | Genetics 



